For Reference 


NOT TO BE TAKEN FROM THIS ROOM 





Gx wpnis 
UNIMASTTARIS 
HARERTAEASIS 





BRUCE PEEL SPECIAL COLLECTIONS LIBRARY 
UNIVERSITY OF ALBERTA LIBRARY 


REQUEST FOR DUPLICATION 


I wish a photocopy of the thesis by 
AS +O ° (author) 


entitled | 





The copy is for the sole purpose of private scholarly or scientific study 
and research. I will not reproduce, sell or distribute the copy I request, 
and I will not copy any substantial part of it in my own work without per- 
mission of the copyright owner. I understand that the Library performs 

the service of copying at my request, and I assume all copyright responsi- 


Digitized by the Internet Archive 
In 2022 with funding from 
University of Alberta Libraries 


https://archive.org/details/Johnston1972 














THE UNIVERSITY OF ALBERTA 


TRAVELLING WAVE AND RESONANT MICROWAVE HEATING APPLICATOR DESIGN 


by 


DONALD ALEXANDER JOHNSTON 





A THESIS 

SUBMITTED TO THE FACULTY OF GRADUATE STUDIES AND RESEARCH 

IN PARTIAL FULFILLMENT OF THE REQUIREMENTS FOR THE DEGREE 
| OF DOCTOR OF PHILOSOPHY 


DEPARTMENT OF ELECTRICAL ENGINEERING 


EDMONTON, ALBERTA 
FALL 1972 








cout ATRSAIA BO: TIE RNIN FHT 


WOI230 AOTADITIVGA OMITASH AVAWOMOIM THANO2IS GMA AVAW OMTIJSVART 


2r2qMt A 

- HOMAD2IA OWA ZaiduT2 JTAUGARA AO YTIVOAA NT OT GAaTTIMAuA 
aaa aT 807 2THAMINTUNAA AKT IO THRMASTAIO SAITHAS AT 
YH902041H9 30 907900 0 


“at ¥ ‘ » 
9.4 le 


‘ 


ac a 


THE UNIVERSITY OF ALBERTA 


FACULTY OF GRADUATE STUDIES AND RESEARCH 


The undersigned certify that they have read, and 
recommend to the Faculty of Graduate Studies and Research, 
for acceptance, a thesis entitled Travelling Wave and 
Resonant Microwave Heating Applicator Design submitted by 
Donaid Alexander Johnston in partial fulfillment of the 


requirements for the degree of Doctor of Philosophy. 


if Gs 


; arash: W ata guT 
| HIRATCIA IM SSPAITS WADARS VO-VT AAT 


bus abner avn vadd teddy yitiveo benpteysbay ett 
~faeseeR Sew esi butd eteubs Yo viluset sd? OQ) Breamozet 


bas sve catifaverT.beltting 2fzed) & ,S0nRFga906 70% 





wi betdivelu2 onbeat wifkeotlaqe Oni sesh syvawet: PM Snsno2sh - 
oid Ye tment lifut fabtves at notaiodl. veinexslA bi saed : 


7 
‘ 
vanocalnta te «ctond Yo 4setneb sf? we? 2tnemeiiwoest 





III 
ABSTRACT 


The dominant problems of travelling wave and resonant 
microwave heating applicators are uneven heating and reflected pow- 
er respectively, although both problems occur with either type. 

It is shown that a travelling wave applicator using the TEM 
rather than the more common TE19 mode gives, in principle, improved 
uniformity of heating for a wide range of web materials. Wave imped- 
ance and attenuation functions for the TEM rectangular structure are 
derived using the method of moments. The results show that a ten-to- 
one range of attenuation control is possible by varying the gaps be- 
tween the sides and the centre conductor. This also permits control 
of reflected power and leakage. Efficiency and uniformity of heating 
in a batch-heated web in a travelling wave applicator are analyzed, 
using the initial equilibrium temperature of the web as an independ- 
ent variable. Vapor pressure, standing wave effects and changes in 
wave velocity are included. It is shown that batch heating theory can 
be used to give bounds of inhomogeneity in moving wet webs. This per- 
mits the choice of an attenuation factor. 

Reactive coupling of energy into resonant structures is ana- 
lyzed. Reflections result when waves in the coupling element and the 
cavity travel in or near synchronism. This type of reflections, as 
well as the heating inhomogeneity, are reduced by a new slow wave coup- 
ler design. A cavity applicator providing multimode operation in a 
small volume and utilizing a system of interconnected TE on resonators 
is designed. Tests with the slow wave coupler indicate that the 


frequency pulling effect can be minimized. 
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Z - normalized Smith chart variable (-) 
a - real propagation constant (nepers/m) 
a - coupling coefficient (-) 
8 - numerical analysis acceleration factor (Ch.3) (-) 
8 - imaginary propagation constant (1/m) 
8 - coupling coefficient (Ch.5) (-) 
y - sign variable, numerical analysis (-) 
r - reflection coefficient (-) 
§ - Dirac delta function (-) 
e - dielectric constant (F/m) 
_ Dielectric constant of free space (F/m) 
ep ze - relative dielectric constant, real part, imaginary part (-) 
Behr " of water (-) 
n - waveguide wave impedance (2) 
Ny 7 free space wave impedance (2) 
6 - temperature hal 


6, - equilibrium temperature (°c) 


(m) 
(m) 
(m2) 
(1) 
(~) 
(m\2aqor) 
he) 
(-) 
{m\T) 
(=} 
(-) 
he) 
(-) 
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Aus latent heat of evaporation of water 

do? dg - wavelength 

Upeh - permeability 

&2 - electric field concentration factor (shape) 


o - water filling factor in web 


o - charge density 
t - Normalized time 
6 - voltage 


> - average voltage for paired charges 


ave 
diff - voltage difference for paired charges 

v2 - electric field concentration factor (reflections) 
w - (radian) frequency 


Wy 7 resonant frequency, n'th mode 


i aes waveguide cutoff frequency 
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CHAPTER ONE 


MICROWAVE POWER APPLICATORS 


lade Introduction 


Efficient, high power magnetrons were developed for World 
War II radar technology. The power carried by a radar beam could 
be felt as heat and it was natural to wonder if this effect could 
not be used. One of the first practical heating applications was 


a microwave oven'°9) 


In the intervening years between the first 
oven and the present, further applications and methods have been 
developed, and a small industry based on the heating effect of 
microwaves has been established. 

In its production, microwave energy is far removed from 
that stored in the original fuel, particularly in terms of the 
capital expenditure involve’. The usual production cycle starts 
with conventional electric power which is rectified, giving a high 
voltage pulsating direct current and then converted to microwave 
power by a device such as a magnetron. Magnetron depreciation, 
at one to five cents per kwh (69) is presently the largest cost 
factor. Although efficiencies are high, over 75% for many power 
sources, some power is lost at each conversion step. For these 
reasons microwave power is expensive and its use must be justified 
by factors such as improved product quality, increased capacity of 
a manufacturing process, or elimination of a difficult heat 
transfer problem. The cost of microwave power is steadily 


decreasing as the technology in the industry improves. 
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The uses of .microwave power are now quite divansesbe). 


For 
example, the food industry has employed it in drying potato chips, 
proofing doughnuts, cooking chickens and roasting nuts. Institutional 
outlets can use microwave ovens for quick reheating of almost any food, 
usually with better flavor than can be produced by any other method. 
The forest products industry has employed microwave power in drying 
wood, veneer, paper and special products such as baseball bats. Other 
uses by the plastics, photographic, auto, drug and other industries 
are listed in books (5!) and review articles(2) > (69) on the subject. 
There is now enough experience with microwave power that it is 
becoming evident that some standard techniques should be refined. 
Industry has concentrated a major part of its development 
effort on power sources. Radar magnetrons were designed to produce 
large but brief pulses of power. The original task was to produce 
devices which generated a waveform carrying more average energy. 
The resulting rugged magnetrons have served wer (12) (16) | but there 
are still two problems associated with these power sources. The 
most important is the cost, and secondly, their sensitivity to 
loading conditions. The sensitivity to loading conditions can be 
improved by using devices such as klystrons and triodes in lieu 
of magnetrons. Efficiencies of 70% are quoted for klystrons, at 


(53),(29)_ 


power levels competitive with magnetrons Triodes do 


not operate well at frequencies comparable to the inverse of the 
electron transit time. Nevertheless they have been advocated, 


at efficiencies of about 25%, to overcome the loading problem'22) , 





wor. aersvrb atinp won s16 sewoq sveworotm to 2ocu OAT a 
caqita ofstoq gate af tt Peete 2k CERO, Sak 9 qmsxa 
Ienobiudttent .2due paidesoy bas enodofrts pntvoos .2tumiquob par toorg 
boot ys azole to entieadey Aotup tot 2rtevo oveworotm sau nso evoltuo 
-boittom veto yns yd beauborg od ma asd rove! retisd Adtw yl Tsuen 
oniyah nt vawoq ovawornim beyolqme asd yyteubnt edouborg Jes70¥ aT, 
yott0 .2d6d [fiedoesd 26 foue etouborg [stooge bas Y9qsq .199N9¥ _boow, 
zetitzubnt veito brs puyb .osus .otdgsipotorg .zottesiq ods yd 202u, 
toatdue ott no (22)«(8) soropsae worver bas ‘steed at beset ove 
at ti tant yewoq eveworstm rttw sonatyaqxe devons won 2t stodT 
-bonttsy ad bfuore zouptatioot bisbast2 amoz tadt Jnsbive pntmoosd 
tnemqofeveb ett to Ivsq yoLsm 5 bosersnaonoo 26d yrt2ubnl vali 
—goubo%q of benpizeb svsw 2qorssnpsm websh .z99° Mee 1SwOg NO srotts 
souho1g. oF 26w Azed Isntpito ody .towoq to 2e2Tug Yotyd Jud. opysl 
LYpVano 9psy9Vs som oniytys. motevew 6 betsvetsp doidw 29otveb. 
oroitt tud . «(SI prow bevrse oved enoxsanpsm bepews entsfuest odT. 
aiT .29o1w0e tewog szart dtiw betstoozes emefdorg owt If tt2 916, 
ot ysivistanee thedt <yfbnoose brs .3203 of} et Inssvogmt Jzom, 
od nso enotttbnos pntbsol ‘od ystvistanee aT .2norstbnos entbsol, 
ugil nt 2aborat brs 2noviayls 26 dove esotveb pnizw yd. bevovgmt: 
ts <2norden la bail on 216 XOX to astonatoftta .enorsenpsm. to. 


ni 
7 ® : 
za 7 


7 
an 





At the present time there is public concern about the 
effect on the environment of still another form of escaping man- 
made energy. Government agencies now control designs for leakage 
of microwave radiation (4) (49) (58) | Work on biological effects 


of uieromaves ‘ce2 


should help in setting exposure standards. This 
concern is timely for the industry because the heaviest usage of 
microwave power is yet to come. 

The communication industry is a far more advanced user of 
microwave technology than the microwave power industry. From 
standard methods in radar and communication engineering, the micro- 
wave power designer can obtain the underlying electromagnetic 


theory and the technology to transport power from place to placese2z 


(14), (56) However, he will find problems which are not encountered 


in communication systems. For example, most antenna techniques 
concentrate on the far f'2lis, while the entire task of heating a 
load is usually done by the near fields, which in the absence of a 
* 

load are mainly inductive and store energy, rather than radiating it. 
Microwave power also involves other disciplines such as thermodynamics 
which are not treated in communications texts. Publications now exist 
which deal with unique aspects of microwave power (37) 2(51) s(54) 

Magnetrons are susceptible to damage from reflected power. 
Inside a magnetron there is a slow wave structure, around which the 
forward travelling wave propagates in the same direction as a 
by? The electromagnetic field of a dipole can be broken into 
components proportional to 1/r", where r is the distance from the 
centre of the device. Far fields may be defined as those for which 
n= 1. These fields carry energy away, while the near fields, for 
which n > 1, contain stored energy. Alternative usage defines the 
far field region as distant enough that variations in r affect the 
relative phase, but not the amplitude, of fields originating from 


different parts of the antenna. In this usage the near fields would 
contain some radiating energy. 
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revolving electron cloud, which acts as an intermediate energy 
carrier. With the normal type of coupling between the slow wave 
structure and the external transmission line, the reflected power 
is divided equally into waves which travel both forward and back- 
ward around the internal structure. The backward travelling wave 
accelerates some of the electrons rather than absorbing their kinetic 
energy as the forward wave does. Back-bombardment of the cathode 
results, overheating it and causing sputtering of its active surface, 
thus shortening tube life. Most devices also fail to oscillate 
at some loading woharevans (7 

To limit the effects of reflected power on magnetrons, 
either the magnetron must be made ruqged or the rest of the system 
must be designed carefully. Magnetron manufacturers have not 
always been able to assume that the VSWR (voltage standing wave 
ratio) will be low. The higher the power level of a magnetron, 
the more difficult it is to design for high VSWR. In present 
microwave ovens the unpredictable VSWR necessitates a sacrifice 


in magnetron efficiency to extend its traiOds 


One low power 

(750 Watt) magnetron has been designed this way to tolerate 80% 

reflected reer! ho At 18 kW a directionally coupled output 

has been used.* Automatic matching devices are available, but 

do not have the response speed to handle sudden faults if the 

magnetron power level is sufficient to cause it to arc iternariy’©°) . 

These attempts have not yet solved the problem satisfactorily. 
Applicator design can degrade product quality as well 


as the magnetron life. In most processes in which microwave 





Py See JMP 4(2) p85 (June 1969). 
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energy is used, the final product should be evenly heated. Standing 
waves are inherent in almost all reflections*, and hot spots occur 
in a load at places where the electric field components are in 

phase and add, while cooler spots occur where they cancel. Another 
problem, skin effect, occurs in loads which absorb energy too 
rapidly, as the heating is concentrated near the surface of such a 
load. Both effects are treated in more detail in electromagnetic 


theory texts (65) (54) 


and microwave power texts » and must be 
controlled in most applications. One of the advantages of micro- 
wave heating is that energy is absorbed by dielectric losses, 
becoming thermal energy inside the load. This is commonly called 
bulk heating and with care, it can be used to produce a rapidly 
and uniformly heated product. 

In travelling wave applicators, heating is performed by 
a wave propagating in a waveguide structure. The load is a thin 
dielectric sheet, or web, which is inserted into the waveguide. 
Part of the energy propagates inside the web, giving bulk heating. 
However, as the energy travels, some of it is lost as heat in the 
load and the wave diminishes in amplitude, giving uneven heating 
as in skin effect. Only part of the energy is carried in the diel- 
ectric and the attenuation depends on the ratio between it and that 
carried in the rest of the applicator. Conventional design of these 
applicators does not allow for much control of attenuation** ; 
some loads may be heated nonuniformly while other lighter loads 


may give low efficiency. Reflected power is not usually severe. 


f An exception could occur if cross-polarization should result. 


Pest: Frinaing field applicators, which heat the load outside the 
wavequide, provide good attenuation control but may be dangerously 
leaky. See Chapter 2. 
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In Chapters 2, 3 and 4, a technique is found for obtaining 
attenuation control. 

In resonant applicators the energy is injected into a cavity 
and stored until the load (or the walls) absorb it. Standing wave 
effects in such an applicator are unavoidable, but there are ways to 
prevent them from producing unacceptable product quality, scorched 
spots, for example. Multimode cavities are commonly used, large 
cavities with many modes which are resonant near the same frequency 
and which have differing field patterns. If energy can be coupled 
into the fields of several modes, rather than only one, the result 


44) ,(38) An additional 


is a more uniform electric field distribution! 
device for obtaining uniformity is the mode stirrer, which rotates or 
oscillates at a rate that is slow compared to the frequency of the 
power, but is fast compared to the heating rate. It can change the 
modal patterns throughout the cavity and as it moves, the electric 
field inhomogeneities are moved pseudorandomly over the load. 
Conventional resonant structures have been designed to improve 
heating uniformity, to prevent leakage of energy(46) | and to provide 
many resonant frequencies in the band in which the magnetrons must 
operate. It is not a simple problem to avoid reflected power when 
dealing with an applicator which is literally full of multiple reflec- 
tions and which carries many different modal patterns. Some work has 
been reported (7) »(23) »(32) (38) (52) (75) iu¢ the basic theoretical treat- 
ment of this coupling problem has not been published. In Chapters 5 and 
6 of this thesis, cavity coupling theory and a working coupling structure 


are described. Both reflection coefficient and heating uniformity are 


improved in an operating system. 
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1.2 Applicator Parameters 


In this thesis the dielectric constant is given by e = 


Ey Key. ~ Jey); where ae is the usual relative dielectric Coretaneey . 


Then Maxwell's equations, in phasor form, have no explicit conduction 
term. Both conduction effects and dipole effects are contained in the 
Poyntina's theorem can be derived by intearating the diverg- 


mh (14). 


ence of the vector ExH over a volume This gives 


Jv Bxti* cumpec Maseedi* Sheds) toe. Besthotubrd aly (1.1) 


where E iatthe RMS electric field vector, D = cE is the electric flux 
vector and the asterisk indicates the complex conjugate of a quantity’. 
Similarly H and B = uH are the magnetic field vector and flux vector 
respectively, where p is the permeability. V is the volume, S is its 
surface, n is a unit vector normal to S, and w is the radian frequency. 
Taking real and imaginary parts, 

Re f EH (-n) dS = equl ef E.E* dV (1.2a) 


Im f ExH-(-n) dS = j2uf (uH:H*_ E-E*e 1) dV 
S 
ee 


The left-hand side of Eq. 1.2a is the power flowing into the volume 


(1.2b) 


and the right-hand side is the dielectric loss, or heating effect, in 
Watts. The right-hand of Eq. 1.2b is the difference between the two 


types of energy stored in the volume, multiplied by 2w. From Eq. ].2b 





gsMany texts also use E and H as the magnitude of E exp(jwt) etc. Here 
they are the (time) RMS values, the magnitudes multiplied by 1/v2, 


with the factor exp(jwt) understood. 
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it can be shown that the energy stored in a lossless cavity is equally 
divided between the electric and magnetic fields. yu is assumed to be 
real and equal to its free space value, which is true for the lossy 
dielectrics usually used as loads. 

The bulk heating effect occurs inside any lossy dielectric 
object, that is, one for which Ep + 0. The power absorbed is given by 
EG 23) ea Tor the whole volume or any part of it. If the applicator 
is a resonant structure and is lossless the integration may be perform- 
ed over the whole structure. The Q of the structure is defined as the 
ratio between the energy stored in it and the energy dissipated per 
radian. If the energy is equally divided between field components and 
Ep is constant, 

E-E* dv ct J OE-E*av 


[NS Sre ee ee ee (1.3a) 
felE-E* dv f E-E*dV+ f[ c 


and if We. << (V - AV), 


e" J E-E* dV 
Geen AV (1.3b) 


iF E-E dV 
where AV is the volume of the dielectric and V is that of the entire 
applicator. V - AV represents the empty volume. 
The three factors which determine the absorbed power are E,9 
the volume of the load and the electric field strength inside the load. 
If En is constant, the latter two factors can be conveniently expressed 


in terms of the filling factor F and the field concentration factor ¢€. 


F is normally defined by 


8 
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Poe AY (1.4) 


st a (1.5) 
AV B<E" dy ; 


This definition makes — equal to the ratio between the RMS (spatially 
averaged) electric field in the load and that averaged over the entire 
structure. | 
This pair of definitions simplifies Q: 
i ey Fee 


if the applicator is lossless*. In Chapter 5 Eq. 1.6 will be extended 


(1.6) 


to include the losses which occur in all practical applicators. 
Travelling wave applicators are waveguide structures loaded 
by dielectric sheets. In this case Q is not a very meaningful quantity. 
What is more important is the ratio between the power absorbed in a 
section of waveguide and that transported through it. From Eq. 1.2a 


the power transported through S by a forward travelling wave is 
= Ly. f.fe 
P “ Ic E.E* dS (ied) 
where n= E/H and is the characteristic impedance of the applicator. 
The power absorbed in an infinitesimal lenath dz of a load with 
cross-sectional area AS is 
aP J 4, gfded. ep 
(5) dz = -WE GE), ( ase :E dS) dz (1.8) 


where AV= AS dz and E is nearly constant in dz. Dividing Eq. 1.8 by 


Eq. 1.7 and using Equations 1.4 and 1.5 we obtain 


ner contains the frequency dependence. 
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are -We e"nFee 


off (1.9a) 


v— 


or ”) 
P = Ps exp(-we ,¢ onFe z) (1.9b) 


where Pe is the incident power at z = 0. The attenuation coefficient 


a (nepers per meter) in 


P = P. exp(-2az) (419c) 
is given by 
a = (1/2)we,etn Fee. (1.10) 


Thus heating uniformity and efficiency can be changed by varying the 


- of the load or by varying 


size and e 

In all previous cases & has been evaluated from electric 
fields existing in the loaded applicator. This presents a difficult 
computation problem if the loading conditions are variable. The field 
patterns of the unloaded applicator are usually easier to obtain and 
perturbation equations may be derived by taking vector products of 
these fields and the fields which would exist after the load is in 
place‘*). If F is small the effect is localized. Only the fields near 
the load will change, and they may be found, by using the local 
boundary conditions. Appendix I gives further information. 

The load in travelling wave applicators is usually a thin 
sheet or web. If F is small and the web edges are parallel to the 
electric fields of the unloaded applicator, the fields will be almost 
unchanged by the load, while the flux D is multiplied by «. The mode 
is called an LSE mode (longitudinal section electric) in this case 
while if the magnetic fields are parallel to the load the mode is 
called LSM (longitudinal section magnetic). LSE modes give higher 


attenuation than LSM modes, where € is reduced to &/e, by the load. 
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In Figure 1.1 two common waveguides are shown. Thin webs 
inserted as shown will result in LSE modes with the same field patterns 
as the dominant modes. Let B be the total length (meters) exposed to 
the microwave energy in a cross-section and let t be the web thickness. 
Then for nearly flat webs S = Bt/2. The field concentration is nearly 
constant over the thickness and as long as the perturbation is small 
(See Appendix I) & is independent of t. Then for any length 


ee inna a shop i RS (1.11) 
eee ree 


S 
This is the form used in the perturbation analyses in the next three 


chapters. 
t 
; B 
Z 


Figure 1.1 Cross-Section of Web in Applicators 
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CHAPTER TWO 
ATTENUATION CONTROL IN TRAVELLING WAVE APPLICATORS 


This and the next two chapters form a section dealing with 
web applicators. In this chapter conventional applicators are examined 
for attenuation control. Some new travelling wave applicators are then 
described and the rectangular TEM structure is singled out as the most 


promising, a detailed analysis of which will be given in Chapter Three. 


Eh Conventional Techniques 


The use of rectangular TE19 waveguide is still very common 
for microwave heating of webs. The web passes through slots in the 
center of the broad side of the waveguide, as shown in Figure 2.1. 
Each time the guide passes across the web, part of the energy is deposited 
in the dielectric. To use the remaining energy, the wave must make 
several passes, as shown in Figure 2.2. This is a crosscurrent appli- 
cator, i.e. the web moves forward slowly while the microwave energy 


(35) 


crosses it (See Chapter 4.) 
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Figure 2.1. Cross-section of TE ) Applicator 
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Usually such heating is used for moisture removal. Bulk 
heating of a dielectric delivers the latent heat of evaporation to 
the water without conventional heat transfer through an insulating 
surface layer. The air supply shown in the figures is used for moisture 
transport; it need not be heated if it is dry and moving rapidly. 
Furthermore, the walls are cold and will condense the water from 
overly moist air. Small turbulent jets impinge upon the web at right 
angles and large holes are provided for easy air escape. Any liquid 
on the walls must be blown out or fall out of the bottom exits. The 
guide must not drip, since this could cause arcs. Further detail is 


provided by Okress (2!) 





Motion 


Figure 2.2 Bottom Half of TE) 9Applicator 
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As the wave crosses the web, it loses energy. Altman 2) 
uses perturbation theory (See Appendix I) to obtain: 

a = we oe, FastE,: E) ds (2.1) 

25, (E,* H)-n ds 

where a is the attenuation coefficient in nepers per meter, Ee is the 
unperturbed electric field vector, E is the perturbed field, and n is 
a unit vector perpendicular to the cross-section. If the web does not 
change the fields Eq. 2.1 gives the same expression as before, 
a = (Se ,euFeon. This is a general description; it shows that reducing 


either F or 52 will improve uniformity at the expense of efficiency. 





Figure 2.3 Currents and Dimensions for TE19 Applicator 


Applying Eq. 2.1 to the standard TEy9 applicator, 
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factor Fy is given by t/a of Figure 2.3. The guide wavelength dg is 
greater than the free space wavelength doe an effect which gives increased 
Og near cutoff. 

In a line along the centre of the broad wall of the guide the 
wall currents flow only in the z direction. At other x locations in the 
broad wall they will cross the slot, producing displacement current, and 
from it, radiation or leakage. This fixes x/a to 1/2 and ss TO. 2. 

There is no acceptable way of controlling E unless other structures 
are introduced into the applicator to change the basic mode tt!) 

The width a can be changed to control F but cutoff results 
if it is less than d/2- It should not exceed 3A ,/2 because other modes 
can be generated in wide guides. Changing b deposits more power in 
less web, leaving Og unchanged. Thus simple TEq9 applicators do not 
permit good attenuation control ¢49): 

Another type of applicator is the fringing field applicator, 
where a field component exterior to the structure intercepts the 


q(19), (42) Very little power is carried in such fields end they 


loa 
die away rapidly with distance from the applicator. By using multiple 
resonator sections, very high fringing field strengths may be 

obtained and since the bulk of the power is not applied to the load, 
the attenuation may be gradual. Fringing field applicators are 
distributed structures, so they will accept a wide variety of loading 
conditions with little reflected power.* When operated unloaded, 

by This effect is discussed in greater detail in Chapter 5 where 
slow wave structures are adapted to cavity coupling. The structures 
used in that chapter are fringing field applicators applying power to 


an intermediate structure, a multimode resonant cavity. Phase control 
is also important to that problem, while it is nearly irrelevant here. 
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power is radiated away and the reflection coefficient can be controlled. 





Figure 2.4 TEM Fringing -Field Applicator 


Some sample structures are shown in Figures 2.4 and 2.5. 
In the first case, a stripline has been formed into meandering 
transmission line and the loag above it absorbs power from the fields 
emanating from the top of the conductor. Before terminating on the 
ground plane, these flux lines are intercepted by the load, and the 
further the load is from the conductor, the fewer flux lines reach 
it. Most of the power is carried below the stripline. In the second 
example, one side has been removed from a meandering TE19 waveguide 
and the fields inside spill out. Air may be blown from below for 


flotation and moisture transport. 





Figure 2.5 TE, 9 Fringing-Field Applicator 
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The fact that fringing fields die away rapidly implies 
that only the surface of a thick material will be heated. Therefore 
fringing field applicators have been considered for drying the paint 
strips on pavement (24). although this has not proved feasable. Radia- 
tion of unused power has the advantage that it can protect the magnetron 
from reflections, but the result can be accidents. There are also hot 
spots in almost all of the fringing field structures, but this problem 
may be overcome by adjusting the dimensions of the applicator*. 

Because of their insensitivity to loading conditions these 
structures are replacing the TE16 applicator in practice. However, 
many of their advantages can be more safely attained in a closed 


applicator, as will be shown in the rest of the chapter. 


owe Circular Coaxial Applicator 





Figure 2.6 Circular Coaxial Applicator 





r Pat. pending in 1971, Genesy's Systems. 
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A simple TEM alternative to TE19 waveguide is the circu- 
lar coaxial structure of Figure 2.6*. Its electric fields obey the 
well known relationship: 


Sh A ae, 


r a. r (2.3) 
r Tog, (F 


where E. is the radial electric field, a, b and t are as in 


Figure 2.6, r is a radial unit vector and r is the distance from 


(56) 


t 
the centre If t is small, 


E. (2.4) 
Substituting in 2.1, we obtain 


el b ’ 2 
coe te) le 2wf t Creeying dr (2.5) 


An b 
2f> 1 
a Tp (~ Atay ) -2nrdr 


Ege pny tw (1 + b°/ab + a*/ab) 
30 a+ D 


(Ee 


or 


e"n we. (1 + b/a + a/b) (2.6a) 


a ofr" 


* 
where the subscript c refers to coaxial lines, and 


cig ael 1h (2.6b) 


©. cnla-t b) 


The filling factor ra is defined as in the previous chapter. If we 
separate out the field concentraction factor, we obtain 

2 

— =.(1/3)(1 + b/a + a/b) (ZY) 


ba Another alternative is the twinline, in Appendix II. 
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This becomes unity as the characteristic impedance becomes very small. 
In all TEM structures the currents are parallel to the 

direction of propagation. Therefore, a longitudinal slot may be 

located anywhere with no problem of radiation from displacement current 


crossing it. 





Figure 2.7 Attenuation of Simple Structures (From 40) 


Comparing the TE rectangular waveguides to 50 ohm coaxial 
line, we obtain 


a re We 
<= 0.623 (F=)(>>) (2.8) 
g g “g 


In Figure 2.7, computed results for attenuation in a 0.24 kg/m? 

web at about 25% water content are shown. It has no cutoff frequency, 
so the TEM structure can give higher attenuation because re can be 
made large by making the structure small. Cutoff effects prevent 

this in the rectangular waveguide, since t and a determine Ege 
Circular coaxial line provides some attenuation control, and the 
centre conductor forms an obvious air conduit. However, it is not 


simple to build. 5? is a constant if the characteristic impedance 
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is kept constant, and it is close to unity. 

A section of 50 ohm, 1 5/8" circular coaxial line with a 
1/8" slot was constructed. Its walls were 0.08" thick. Using strips 
of very wet paper, the attenuation factor in Figure 2.7 was confirmed. 
For comparison a TEq9 WR284 waveguide section with the same wal] 
thickness and slot width was prepared. This too gave the predicted 
attenuation. 

If a wet web is used, energy propagates out of the 


(69) To check the leakage of the 


applicator in a dielectric wave 
two applicators, they were loaded very heavily with wet paper. A 
laboratory generator applied 100 mW of forward power to each 

applicator. With a Narda 8110 leakage detector it was found that 


5) This means 


leakage was proportional to the ratio of Ot and a 
that leakage depends on the electric field intensity in the web 
region, and not on the mode. When metal dipoles were used, which 
is an unusual practice, the coaxial line gave significantly more 
leakage, since the dipoles were capacitively coupled to the centre 
conductor. This effect was not observed for even the wettest paper. 
Thus for reflectionless dielectric loading, the escaping power is 
proportional to (<*P_)/S where Pe is the forward travelling power. 

A TEM applicator has been found which does what a TE, 


one will. It also permits control of PS but — is nearly constant. 


Zoo Variable Cross-Section TEM Applicators 


In addition to control of F, by scaling the dimensions, 
we may control <2 by placing the web in a region where electric 
fields are Paduced to? « To do this, non-circular cross-sections 


must be used. Because the currents are not transverse, virtually 
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any structure comprised of two or more insulated pieces of metal 
may be slotted longitudinally and used as an applicator, provided 
that it bears no higher order modes. Many possible systems provide 
the reduced a2) giving reduced attenuation, leakage, and even 
reflections, since low — implies that most of the energy does not 


impinge upon the edge of the web. 
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Figure 2.8 Elliptical TEM Applicator 





The simple circular structure can be distorted and made 
elliptical as in Figure 2.8. The region of high power transport 
is the narrow gap between the major and minor axes of the inner 
and outer conductors respectively. The web passes through the low 
field region between the other axes of the conductors. However, 
Be is increased in the process and a is not reduced as rapidly as ee 
A much easier structure to construct is shown in Figure 2.9. 
A rectangular centre conductor * is used with an elongated rectangular 


* In the rest of this thesis it is assumed to be square, but 
if its shape is adjusted, added control of &“ is possible. 
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Figure 2.9 Rectangular TEM Applicator 





é is the same as in the elliptical 


outer conductor. The effect on & 
applicator. The corners must, of course, be rounded to prevent 
arcing. 

A first order analysis may be made by treating the upper 
and side gaps as parallel plate capacitors and ignoring curvature of 
the fields. This is valid for any frequency, since there is no cutoff 
frequency and therefore no dispersion in TEM structures. If we assume 
that w>>c, we can assume that the voltage gives the electric fields 


directly, or o = (E,)d = (E)c. For small perturbations and square 


centre conductors, we may expand 2.1 as follows: 


ill bars rs 
_ Ag®Egey iF (7) tdy 
Se ern ead Tt a 
C 1V ry 
2 vy “wdx + 2 Vv) wdy + E|ds 
byl? fe (q) wdy ejetel 
If d>>c the first of the lower three integrals dominates and the 
power is carried mainly above and below the centre conductor. The 
attenuation is then proportional to c/d. 
This result is for the asyptotic case: low characteristic 


impedance and very small c/d. When the gaps c and d become large the 


electric field is no longer constant and the upper integral becomes 
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a constant because the fields diverge at the side of the centre 
conductor. It is expected that the resulting minimum a may be 
reduced by making the centre conductor higher than it is wide, or 
by making the face of the centre conductor concave. 

The leakage of such a structure is dependent on S while 


attenuation depends on both ee 


and F. Leakage is then proportional 
to (c/d)?, and this applicator can give reduced leakage while 
providing the same attenuation as a given circular coaxial or 
rectangular TEq9 applicator. 

If the characteristic impedance of an applicator undergoes 
sharp changes, as when a heavy load passes through a high field region 
or a is changed, these discontinuities will reflect power toward the 
magnetron. If there is an absorbing web between the reflection and 
the power source, the problem is reduced, but this cannot be guaranteed. 
The characteristic impedance change is proportional to €& , so reduction 
of —€ also controls the problem of reflections from the web edges. 

Similar results have been demonstrated in waveguides 2) » (41) 
by reducing the fields in the centre of the TE, rectangular structure 
by providing two constricted regions on either side of centre. The 
power was conveyed in these regions, reducing ge by the same principle 
as in the TEM structures of Figure 2.8 and Figure 2.9. The biggest 
difficulty in either mode will be arcing, but with rectangular TEM 


structures arcing will be easier to control since wider regions are 


provided for power transport*. 





= A simple TE19 waveguide can carry more power than either 
structure. See the example at the end of Chapter 3. 
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CHAPTER THREE 
NUMERICAL ANALYSIS OF THE RECTANGULAR TEM APPLICATOR 


In the previous chapter it was shown that control of the 
parameter ef is possible in TEM structures. Approximate behavior 
of attenuation for the rectangular version was derived from 
simplifying assumptions. In this chapter Ee will be evaluated more 
accurately and the previous expression will be verified. The 
perturbation technique will be retained, allowing the fields to 
be calculated for the empty applicator. The loaded applicator is 
assumed to have the same fields, which is approximately true for 
the electric field if the web is thin and is parallel to the 
electric fields calculated for the empty applicator. 

The assumption of parallel plate capacitor behavior wil] 
be replaced by the method of moments. In this method it is assumed 
that the charges existing on the conductors can be accurately 
represented by a screen of infinitely long line charges. When these 
charges are obtained, it will be possible to present the designer 
with charts, from which he can select rectangular TEM shapes to give 


desired E2, F and characteristic impedance AGIA oe 


xa Characteristics of the Rectangular Structure 


re can either be derived from the area integration in 
Chapter 1, or from a volume integration over some length of applicator. 
Let AS be the cross-sectional area of the empty structure which will 
later be occupied by the web. Expressed in terms of volumes, if the 


fields in AS are constant over the thickness of the web, 
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(3.1) 


where AV is the volume z'AS. AS can only be changed in magnitude 
by changes in thickness t'. V is the volume z'S and z' is the 
active length of applicator, measured along the propagation path 
of the wave. Similarly Ee can be treated as a power transport 
parameter. The power carried in the total cross-section S at 


unity VSWR is given by 
P. = f(E x H)ds= (4) f(E - Ex) ds (3.2) 
aner ns 


from Poynting's theorem. The power carried in the web region AS 


of the empty structure is given by 


] a a 
Pay mL aete - E*)dS. (Jest 
Then 
2 
e“ = gim (S/AS) (Pe, /Pe) (3.4) 
im fw P+ 


For LSE modes, if the web is found to be a minor perturbation by a 
check such as that given in Appendix I and if the fields are nearly 
constant in its vicinity, the total mean-squared electric fields in 
it can be found by multiplying ae by PeF/n. If pt is known, the 
heating rate and attenuation constant can be found from the power 
level and web properties. 

The electric fields shown in Figure 3.1 may be calculated 
by the method of relaxation (18) ,(27) (47) (60) 4), the method of 
doments (22) It was decided to use an extension of Harrington's moment 


method for economy of computer time as well as for accuracy. The 


relaxation method is iterative, involving a very large number of 


(1.8) 


eas | Wren 

‘sbustnpsm nt bepnsda ad ino hea is RTC 
ait 2t 's bre 2's smlov oat at ¥, 3 seamioid ak zepnats vd 
fidsq notfspaqorq sdt pnols bsavesom KR epeny to dipnaf avtsos 
Syogens1d yswog 6 26 betsevt od m8 ‘3 uivsttmt2 -ovew oid to 
36 2 nofsos2-22019 fstot./ed3 ni betwiso Yewoq edT. .vsdems sq 
yd nowtp 27 WWeY ytiow 


(S.€) | 2b (#3. a) oor apes. A 


Cac 


24 norpey dow ot at bstv1s9 yowoq onT aasiitt 2 ‘tame’. aa 
vd: nevie ef arussounte yiqme ont to 


(€.£) | .2(*4 + Dy bays 


(B.£) >. Gag) (eae) tig = Ss a 


Cert 


6 Yd noitsdwiysq tonim 6 od ot bnuot 2t dow add tt -20b0m ws 


ve A 


ylyean sv ebfott sit tt bas I . eee nt bea ied ar 








26 


calculations. In addition, it results in a grid of voltages 

whose gradient must be evaluated and squared. Numerical 
differentiation can produce very large errors. The method of 

moments produces a set of charges, from which all the other parameters 


may be obtained directly. 





Figure 3.1 Rectangular TEM Applicator 


For matching the applicator to the power source as well 
as for minimizing internal standing waves, the characteristic 
impedance must be known. This is given by Lo = 1/cC where c is the 
velocity of the wave which in empty TEM structures travels at the 
speed of light. C is the capacitance per meter, given by 

C = (la,| + lal )/2 o> (3.5) 
t and qb being the charge per meter on the inner and outer conductors 
respectively. $ is the voltage difference between conductors, and 


the average voltage of the system is set to give equal {q, | and |q,|- 


cpire A Mathematical Technique for TEM Structures 
(a) The Moment Method 


The method of moments can be used for electrostatic systems. 


TEM applicators have the same electrical properties at any frequency 
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if they are reflectionless. Thus we san start from the fact that 
the voltages on the conductors are constant on the metal surfaces 
and find the zero frequency charge distribution which generates 
them. 

To avoid solving integral equations, which involve an 
infinite number of charges, the moment method approximates the 
charge distribution by elementary charged areas. These elements 
have an assumed internal charge distribution or weighting function, 
which is fixed in form but multiplied by a scalar. The scalars have 
the dimensions of charge and convert the relationship into a matrix 
equation for which they form the vector of independent variables. 

In fact, though, it is the voltages which are known. This means 
that a matrix inversion is needed to change this simple relationship 
between charges and voltages into One in which voltage is the 
independent variable. 

The method of moments is an approximation. It is based on 
the fact that the voltage at a distant location* is sensitive not to 
the arrangement of charge in an element, but to the magnitude. 

The approximation can be improved by increasing the number of elements, 
and by selecting the weighting function to better approximate the 
actual distribution or to produce errors which cancel. 

Let a charge density o Coulombs per square meter be 
distributed over the conductors of Figure 3.1. Assume that it 
generates a yoltage d onit 1 volt on the inner conductor. These 


voltages, at any desired location, (x,y,z) can be obtained from o 





be Further from the nearest charge than the width of the 
elementary areas. 
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by electrostatics: 


o(xsy) = fffretey' sz") dx dy dz (3.6) 


Aner 

where r = v(x - x") +(y - y')é + (z- 2')e and the dummy 
coordinates (x', y', z') are the location of the infinitesimal 
charge element of the integration. However, the charge distribution 
o must be found. To do this, we use line charge weighting functions 
(Dirac delta functions), 6(x - x')s(y - y'). The charge sheets o on 
the conductors are then approximated by a screen of line charges 
q; S(x - x')é(y - y'), infinitely long in the z direction. Then 
in matrix tensor form, Equation 3.6 may be approximated by 

oO ssyy) 2p sP ayy") (3.7a) 
and in inverted form, 
gi" pey'y) = Li 
where primes indicate a charge location, and where repeated indices 


q (x' sy") (3.7b) 


i 


indicate summation . Lj is given by 36 ;/2q,. 


Be ek j = Xs and y j =y 


evaluate the self-potential of a line charge, a discontinuity 


;> or in other words if we try to 


results. To avoid this, a strip charge is assumed in evaluating Lis. 
The weighting functions must be changed in this case, and we may use 


strips of width 5 baa 





centered at 


! i] ' ya 2 
(x 52177 j + (y j+] ag cpa 


(x' joy"). 


It is now possible to evaluate the Lis directly from 


electrostatic theory and invert the matrix by digital computer*. 


bad There are inversion routines in the IBM standard routines 
based on Gauss Jordan reduction. One of these was used. However, 
it has been shown that the moment method can give rise to ill- 
conditioned matrices, for which inversion gives escalation of round- 
off errors. This will be discussed later. 
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Figure 3.2 Line Charges 


Let p be a measurement point, a distance r away froma 


line charge of q Coulombs per meter. Then the voltage at p is 
bo = oe a05. 


Ane " 


where r' = Wp? +; 7 . For a single charae the voltage is 


i a (z + /2° + r°) 


oy = gim ee (3.8) 


u 
0 
The infinity may be removed by evaluating the potential for 


equal and opposite pairs as in Figure 3.2*. Then, 


il 
Oo 
Lo | 
+ 
N 
ine) 
8 


ep 





a ztstwicrons af |° 
q: r 
= a an G (3.9) 


In 3.7, if p represents the location of another charge, 


Lod must be expressed in terms of charge j and its image charge, i.e. 


*It may also be removed by using a reference point at a finite 


distance from the applicator. 
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L al 
Ds ab; H Of Sh ane, cofriuc (3.10) 
eReaee, 
The errors in using line charges are evaluated in Appendix III. 


The electric field near a line charge is given by 


eyes (3.11) 
x 
z 
es 
Figure 3.3 Strip Charge 


To evaluate the voltage in the centre of a strip where 
i = j, we use the geometry of Figure 3.3. At the centre, Eq. 3.1 


gives 


dz. (3712) 
Making appropriate substitutions and integrating by parts, we obtain 


b,4 (wb) = 52 [sinn”! (22) + (2) Sinh” (5%) ] (3.13) 
where w is the width of the strip and b is its extent in the +z 
direction. This gives the self-voltage of one of the strips in the 
applicator. Here too, the voltage increases without limit as b 


approaches infinity,* and, a strip of equal and opposite charge must 





* When 2b = w we obtain Harrington's solution for a square area. 
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always be present on the other conductor. 

We may approximate the image charge as we have done for 
the other cases where i # j. A line charge of magnitude q = -ow 
is assumed to be located at its centre. If the matching charge is 


located a distance r away, its voltage for a lencth 2bais, )from 3.8, 


oi ntl-i - oN an (b + A? + r*) sie | (3.14) 


The self voltage of this system of charges is 


dhe thorge Epabrew) (3.15) 
where 
Fi (r,w) = aim pm sinh” | (2) + (2) Sinh CF 


£2 p.(4—— ta Ty (3.16) 
rete 
Taking limits, 2im gen mes 8 2 = gn (5p) gim x Sinh”! ls Bs 
b->co bet yp or Xoo 
and gim Sinh7! (x) = en x + &n 2. 
Xo 
This gives 
F,(rsw) = zim fen (22) + en (5) +1 + an 2} 
boo 
= gn (“) + 1+ en 2 (selza 
and 
> = a (2n at 1.693147). (3.18) 


Eq. 3.18 gives the self voltage at the centre of a strip, and must 


of course be used again with a different w when obtaining the self 


voltage of the image charge. 
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(b) Use of Image Charges 


For the rectangular TEM applicator, each charge occurs 
in the presence of an equal and opposite charge on the other 
conductor. In addition, three more identical charge pairs are 
located symmetrically as in Figure 3.4. Assume that (ow Joo 
is at (Yn 2%) while (ow) 0] is at (y,>-X,) » (o Wi is at 
(- m? Xm) and (ow, 1 is at (-y,»-x,)- The matching charges 


ok 
index of the charge considered, and %,k are O for the positive half 


are (o where n is the total number of charges, m is the 


n+1-m 


plane and 1 for the negative half plane. (O44 om ho ie -O nM, 
: Q k 

and is located at Ly; ie hy, Pe ar geek) ), for a 

total of four images. Then if i # j using Eq. 3.5, a =m, 


b = n+ 1 -m and (3215) 
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Figure 3.4 Image Charges 
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When 7 = j, the self yoltage of this pair must be added 
to those of the image pairs. Using Equations 3.5 and 3.8 


yee, teeter 
eg] (YE-Yngq_j) * G-%p4q_4) 
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Fa 
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“TT Set FOS bay Seer” 
Ay sYyera) + OG 41a) 





tyiab Seererietnsted ir the 7 ine! eb (3.20) 


/ 2 2 r 
(Yi*Yper-a) + OG era) 


It is now possible to invert the matrix Li5> providing that 
the proper distribution of the opposite images is known. This must 
be determined iteratively. As a first approximation, we can assign 
images in a location close to the correct one by simply distributing 
them equidistant from each other. For a magnitude of o.W. in 3.7a 
we may use that of either the charge i or its image nt+1-i, whichever 
is closer. The accuracy of this is improved by the fact that the log 


function is quite insensitive to errors in faraway charges. We then 


let 
L = ] {2 (_tontl-Sy } if r <V7 L =f 
ier alice vee ij —‘isntl-j, “ij 
LL eh ae & or alternatively we may average the two charges, 


ij i,ntl-j 
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r. 
letting Lis = 5 (xe) { en a aot + ...} for each. Both 
Ts 


methods give accurate capacitance results, but to reposition the 


charges later it is desirable simply to use the nearest charge. 
(c) Calculation of the Average Voltage. 


If the charges are not paired properly, the capacitance 
and electric fields will be slightly incorrect. This error will 
arise in practice because the average voltage is not usually zero 
when the charges are paired, but shielding results in some average 
voltage other than zero for the structure. An infinitesimal 
unbalance in these infinitely long charges is sufficient to give a 
finite average voltage to the structure, as is shown by Eq. 3.8. 
This possibility has been eliminated in the final expressions for 
Liz Eq.3.19 and Eq. 3.20. Thus the error has no useful physical 
significance, since it results from violating the charge pairing 
assumption. 

Let us assume a 1 volt potential difference between 
conductors and evaluate the resulting average voltage. To obtain 
the capacitance, let 
qi Lis ay (3.21) 

1 ] 


— - = <= 2 aed 
j bere 5 volt for the outer conductor and 5 > 5 volt 


Then ave 


for the inner conductor, where ove is the average voltage to be 
calculated for the structure. Then 


49 © In+1-4 


and 


; qi |. (3.22) 
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Since 3.21 is a system of linear equations, 


-1 -] 
eel. + 
oS is Fe *o, bave bag Yy (3,23) 
where and u. ; ae) nepal 
ie u are vectors $0. 5 for i > n/2 and ae 5 
J J J 
roGee nie. us = ]. for all jw, Let 
%. = La oe (3.24) 
i ey 
Then 
2 -] 
a 40, * bave Lij me, (3.25) 
n 
We know that ) q, = 0 from charge pairing. Then 
i=] 
: n 
0 = q. + ¢ } 71 
ie a ek 
and we may solve for ave: 
(Cena ey AEE, us) 
b =F = gq Rete u.). 3.26) 
ave et ae et mr ee 
The capacitance is then given directly by 
n/2 -] 
uma (9, ag PR Li u;) [3ce7) 


from which Lo follows. 

These facts can be shown to be reasonable by examining the 
individual elements and their interactions. As expected, it is 
also possible to show that the average voltage is proportional to 


the potential difference. If a ee the average voltage for 
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a matched pair, neglecting other charges in the vicinity, is 


Les Visntl-i Matl-i,i 
bss = -(~—) {2n (—2——) - an (———}} .__ (3.28 
aVeun ck Wi Wat] -7 


The potential difference is 


lw 
1l-i Matt TS 
fee ak sabes Fir pa chalk petty yp eal bea 

iiggee 2m W. wats} )} 





5 ay it aed A | Wer seit oe (3.29) 
Solving for (q./2ne), 
fe fy (3.30) 
“di ff. "i 
>: - of 
‘ave 2 en( ris nv} peniaeny + Qn pfelpted | +2 + 22n 2 
Wi Wat] -i 


All quantities except ints and dare are fixed by the given 
geometry, and unless ener fave #0. The voltages of this pair on 
elements j and n+1l-j in the vicinity are 


q. 
-{ gn (SER ae ae ai) 
ae 
Vig n+1-i 


547 


and the average voltage on pair j, nt+l-j due to i, n+1-i is 


q. Vy saes ta (3,32) 
bave = & GE) Can (SL) ~ ant ST} 
ji "jentl-j nt+1-j,nt1-j 


Again, if the pair of positive charges is spaced differently to the 
pair of negative charges, this contribution to the average voltage is 
not zero, and it is proportional to the potential difference ete oe 
If the charge distribution were chosen correctly, the 
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resultant average voltages would all be the same*. For correctly 
spaced charges, 
n/2 


7 2 
>. = > OY > 3 (3.33a) 
Tave jal 2¥8y4 2 ih ave 


and these averages would all be equal to that given by Eq. 3.26. 
The average voltage in 3.26 for incorrectly placed charges is 


n (3.33b) 


It was found when a computer program was developed both to 
locate and to solve for the charges, that the average voltage for a 
pair of parallel charge arrays was zero. This is predicted by 
Equations 3.28 to 3.32, since the spacing of the positive strips was 
the same as for the negative strips. When a closed coaxial system 
was used, the outer conductor shielded the inner one and the 
applicator had an average potential very nearly equal to half the 


potential difference. 


(d) Correction of Charge Placement 


If the charges are not correctly placed, some error 
remains in the values computed for C and ». Experimental results 
showed the capacitance error to be about lp for simple evenly spaced 
charge distributions on square 47 ohm coaxial line. However, 
after equalizing the magnitudes of the two charges on each pair 
(which makes the situation worse), errors in voltage of up to 20% 


remained. Since the electric field is the gradient of voltage 





Fs For spacings chosen incorrectly, the potential difference 
is still the same for each pair, if the method of moments is used to 
obtain q; as in Equations 3.21 to 3.26. 
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and since the potential difference between image charges is correct, 
its error is less. However, the percentage error in & is double 
that in the fields themselves, and a is also sensitive to local 
charge errors. 

To remedy the error in charge placement, we must examine 
the physical meaning of the line Charges. The ok matrix has been 
set up so that "qj is assumed to be paired by a charge “4. at 


location n+1-j. We then set L. -=0 if i and nt+1-j are on the 


ont] -j 
same conductor. The line charges computed represent a strip with 
magnitude q. equal to o5W,- The location of the image charge is 
nearly, but not exactly, correct and the log dependence of the 

Lay makes them insensitive to shifting of remote charges. Thus one 
may reduce the width of qj to compensate for deviations in q from 


1 
aq 


on the centre conductor fixed, since the corners affect the field 


- Wnt1-j): It has been found practical to leave the distribution 


strengths strongly. One may then make double the indicated 
correction in widths of the strips on the outside conductor. 


Letting qj represent the average charge for the pair qj? 


oa 
9; ey (q, > Wnt1-j)° (3.34) 
Then the new width W. is given by 
t i] ' 
. = eae ; a Spas a) 
w= (1+ B(q, q,;)/45) W (33>) 


where 8 is a multiplication factor; which is different from 2 because 
ave will change slightly. The pairing changes as W. is adjusted 
and the total width must be equal to the distance around the outer 


boundary. 
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When each iteration, consisting of calculation of We followed 
' 
by calculation of q4: has been performed, the matrix Li, may be used 
with q4 to give the voltages on the conductors. In other words , 


] 1 
$- =>5 (L..qi + L. .q :) (3.36) 
Tye 2 ‘is fea 1,nt+1-j ‘nt+1-j 


) may be used to check for i11-condi- 


if 
— 
a 


The error e, ‘ - ¢ 
i ave 


ave 
tioned matrices and/or inaccurate charge distribution. The voltage 
difference is always correct if the matrix is well-conditioned*, and 
this may be used to pin-point the source of the error. When the 
charges are used to calculate the fields, the tangential electric 
field will be zero. Since this implies that adjacent strips are at 
the same potential, the two facts are equivalent. If the charges 

are not equal and opposite, the product q5 = o5W, must be changed. 
This was the criterion used in the program, since 2 is sensitive 

to local errors. 

To use the program to calculate its own 8g, we first 
determine a charge error for each strip pair. Let the error be given 
by 
e, = (q, - q.)/q, (3.37) 


‘e See for example, McCalla, Introduction to Numerical | . 
Methods and al gt pros rameltnd j J. Wiley 1967, p.176. Bostian's 
word of warning\!9) is even more valid for infinitely long ne 
charges than for point charges. Serious errors arose from 7 
conditioned matrices using single precision arithmetic. able 
the voltage difference check showed that typical accuracies : 

six figures could be obtained by using double precision arith- 
metic on the IBM 360/67 system during matrix inversion. pes 
iterative procedures such as the Gauss Siedel method were no 


necessary. 
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where q; is the charge calculated during an iteration, including 
the effect of the average voltage as in the previous section. q. 
is this charge averaged with its image charge. The optimum 


acceleration for an iteration occurs when 


ren iy. 
e = ) ery = 0 (3.38) 


1 ife; was previously positive 
where y; = 


0 if e; was previously negative (3.39) 


The 8 satisfying this criterion would give all e, = 0 if the 


aL. . 


a were constant with respect to q; and Ws. To correct g from 
; 


By to Bo in two successive iterations, we use B to correct q; 


to q. Then we examine the magnitude of the correction. If we 


PZ 
use B, we obtain an error e. . Summing 
2 
> n/2 
Sp Stebride, wi se) 
i=] ] 
and 
hn ne 
e, = a 
2 48, 12 2 


The actual correction for each strip is equal to (e,¥} C079.) and 
1 
if this is proportional to By» we have that 


(ey; zs Ono )/ By = const. (3.41) 
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The desired correction giyes no averaged error, or 


(ery, )/B5 = const. (3.42) 


Then by summing over 7 and dividing Eq. 3.41 by Eq. 3.42, 


n/2 
‘ae 
By = By bt Pia! (3.43) 
n/2 
) | (e, -e, )| 


Eq. 3.43 gives a new acceleration factor for the next 
iteration, or alternatively for correction of the previous adjust- 
ment of w.. If the error e, is less than e) and g is nearly 
constant, convergence is most rapid when By is used in iteration 
ptl. As expected, g is about 2, typically 1.7. When there were 
large errors, the corrections were scaled down until the largest 
adjustment to a strip was about 66.7%. 

Using the methods of this section, small errors converged 
rapidly, i.e. a voltage error of 6% converged to 0.3% in two 
iterations, starting with a g of 2 and recalculating it once. 
Simultaneously, charge errors of 32% converged to 1.8%. The 
method loses speed of convergence rapidly for errors over 50% 
and for applicator cross-sections where the width was very much 
greater than the height. In addition to limiting the corrections 
it was necessary to start from a nearly square structure and 
progressively make the structures wider. Increments of 20% 
were used, which was also a good interval for plotting the 


results, and each corrected distribution was used for the 
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initial strips on the outside conductor of the next cross-section. 

Theoretically, Newton's method could be adapted in matrix 
form, using first derivatives of either >, or tangential electric 
field. When this was tried, it was found that very ill-conditioned 
derivative matrices resulted and the errors were drastic. The 
method based on adjusting total strip charge was much more 


successful. 


(e) Calculation of Applicator Parameters 


When the charges have been calculated and the error 
criterion met, we may calculate 2 TYOMMEGUALIONS rls o.c5 OL0, 
and 3.11 by placing a grid of points in the interior. The total 
fet from all charge pairs is obtained for each point and multi- 
plied by the cross-sectional area around the point, and a summation 
is used instead of the integrals in Equations 3.1] and 3.2. For 


a true TEM structure, 





d % o.W; Z 2 
Smeg ay eh p +43 (3.44) 
Yr r (no summation) 


where x = (X_ - x)» y=(y - y,) and r = x + y as shown in 


p p 
Figure 3.4. The total squared field at point Ps is then given by 


ze 
a; (") 43) te 


Ll n 
(E+E), = (pe) CD 
(3.45) 

If k points are evenly spread over the main area S and 2 points 


over the web area,AS and Aa, is the area around point i, 
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The power carried in the guide for a voltage difference 
daaet is then 


k 
Pai (yy 


(Eu Ey ha. } (3.47) 
(e) j=l J 


Alternatively, 


; (3.48) 


Pee (busee)” pats 
which can be compared to the results of 3.47 to check the computer 
program. To check a particular design for breakdown problems, we 
may use the square root of Eq. 3.45, scaled to some forward power, 
i.e. 1.0 kW using Eq. 3.47. The maximum E in the interior will be 
somewhat less than that at the corners if they are not rounded. This 
means that the value given by a grid is an approximate value only, 


applying to rounding radii of the order of the grid and inter- 


charge spacing. 
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3.0 Results 


Two computer programs were used to produce design maps for 
the rectangular TEM structure. The first was an implementation of 
the mathematics of the previous five sections. This main program 
accepted as input parameters the dimensions of the structure and 
several control commands, and was run for a wide variety of appli- 
cator cross-sections, producing a deck of output cards carrying Ly 
and 5? data. A second program was written to interpolate linearly 
between the members of this set of cross-sections. This program 
was capable of rearranging the data for presentation on various sets 
of axes. 

The main program also printed sheets of data, giving the 
applicator dimensions, Zo? the number of sampling points above and 
to the side of the centre conductor, the error in Poynting's vector 
(obtained by comparing Equations 3.47 and 3.48), the factor 52 for a 
web above as well as to the side of the centre conductor, and the 
filling factor. This information could be used to calculate attenu- 
ation and to predict its accuracy. In addition, the maximum electric 
field was printed for a forward power of 1 kW. This was simply the 
largest field found on the grid of measurement points. It is about 
the same as for a structure in which the corners are rounded with a 
radius equal to the spacing of the charges on the centre conductor 
if the grid of measurement points has the same spacing. The spacing 
could be checked by referring to the printed output, and adjusted 
accordingly. To obtain better accuracy, the charges may be placed 


closer together on a quarter circle at the corners. 
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It is likely that the present method is accurate to 
about 50% when calculating breakdown power levels. 

For the interpolation step, a matrix of applicator 
cross-sections was obtained by varying the upper and lower gaps 
by 20% steps. Including images, 80 charges were assumed per 
cross-section. The inner conductor was square and was assumed 
to be 1 meter on a side, but the dimensions may be taken to be 
any convenient units. The spacing is 1/20 of the inner conductor, 
so fringing fields to the side of it are fairly accurately 
determined. In the computer program, EQ was SEL arOai ern TNs 
means that the charges were not expressed in Coulombs but may be 
converted by dividing by 18 x 10”. 

Figure 3.5 shows the charges distributed over one 
quadrant of a typical applicator and the voltages induced in the 
interior at the sampling points. Figure 3.6 shows the same 
applicator with the magnitudes of the electric fields that the 
main program calculated. The charges on the outer conductor have 
been iteratively relocated and the outer corners have very little 
charge. The corners of the inner conductor have about twice the 
charge density that the centres have. The measurement mesh 
chosen was too coarse to determine the maximum electric field, but 
the corner has a high gradient. Also apparent is the field 
concentration near the centre conductor, indicating that the 
faces of the conductors are not acting like a parallel plate 
capacitor, as they do for smaller gaps. If the applicator is 


scaled in size, Ee, Ly and the voltages are unchanged. The 
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electric fields are inversely proportional to the scaling factor 
and directly proportional to the applied voltage. For convenience 
in computation, the dimensions shown are doubled, while the fields 
are calculated for a normal-sized structure, including all four 
quadrants. 

Figures 3.7 to 3.12 are maps produced by the interpolation 
program. The first four charts are contours of various parameters 
plotted against the dimensions of the structure. With the exception 
of Figure 3.10 for Zo? these maps are not in a convenient form for 
design but they show how the various parameters change as the 
designer changes the vertical and horizontal gaps, corresponding 
to c, and d respectively of Figure 2.9. The third gap, through 
which the web passes, is assumed to be very small. 

A deck of 136 input cards was converted to a set of 6 
matrices of size 16 by 16 whose elements were input or calculated 
data from the main program. The interpolation program searched for 
contours between these points. Each input card was accompanied by 
a card bearing the error calculated in the main program which was 
used to correct —&. Since most of the power is carried in the narrow 
gap, it was assumed that the same percentage error occurs in z2 as 
in Poynting's vector. For almost all the cross-sections the 
uncorrected error was less than 10% and only for very small gaps, 
or Zo of about 12 ohms, was the error even visible on the graphs. 

Figure 3.7 shows the filling factor. All the other para- 
meters have been normalized to the dimensions of the centre conductor, 


so for F, the web is assumed to be as thick as the centre conductor 


8h 


oso6T ihe ait of [snoftiogorg. ents 916 abistt oivsosts 
soretnevies 103 .spsd Tov bona aft of, {snotsogorg ystzert bns 
abfett oft oftiw .bolduob 976 morte enotenamtb aft Hof saduqnes at 
duot its ontbufont Ut UZ boste-femvon 6 Yo? besstuatea 916 
Pa. 

nottsfoquetat sit yd beosubovqg eqsm ove ST.€ oF V.E eo upt? 
avatemeisq evotréy to e1wotnos sts 2tveio Yuet Jevit sxT <mangovd 
notsqsaxs ont astW .swutouyte sAt.to snotenemth sit tenteps bettolq 
+o¥ myo? fnstnevios 6 NT ton 9x5 eqem sess ty sot OF.€ swwptd to 
ott 26 spnsd> 2vetomevsq evotyev sit word wode yedt tud niptesb 
pntbnogzsyvoo .2qsp Istnositod bas Teoftyev add edensdo tonptesb 
dpuords .qse brtdd oT .@.S oyupt? to yfevitosqesy b bas . of 
.ffeme yiov ed od bombees 2f ,2e9e26q dew add dotdw 

& to te2 6 od bedtevAod 26w 2by69 Suqnt act to josb A 
badetual'sp yo tuant ssw etnemsio saonw Of yd Of oste to egotism 


Yot beravsse msip071G nottsfoqrvetat siT .meypoyq atsm ads movt sieb 
xd bsinsqmo276 26w bi62 tuqnt fond .2intog sess nsswted axuosinos 
| ile dotdw mexpoig nism oft nt beteluolso vos oft pnivsed bis 6 

worven oft nt botytso° 2t yswoq sd¢ to ‘teom sonte 2 eialcab i beau — 
ae Nk 214990 O19 ‘omens ompe oft Sond. bamses om 3m ig 









VERTICAL GAP, c 


4g 


10! 






10% 2 3 Cit Tae Ser 2 
10° 10° 10! 


HORIZONTAL GAP, d 
Figure 3.7 FILLING FACTOR NORMALIZED TO CENTRE CONDUCTOR 


Figure 3.7 FILLING FAC; On NU 


- The normalized filling factor FF) is used. AS is the cross- 
hatched area to the side. 
- w taken to be 1 unit, c anc d follow. 





50 


is high. This does not imply that a large gap occurs in the 

centre conductor, but that the computed values of F and Fee must 

be multiplied by t/w where t 7s the web thickness and w is the centre 
conductor height. Thus the plots show the variable F(7). These 
contours tend to become constant for very wide webs. The scatter 

is due to the interpolation program. 

In Figure 3.8 the quantity ge is plotted. For very narrow 
applicators the web is in the high field region. The interior area 
can be very large if the vertical gap c is large. Thus ag can 
increase without limit. The contours tend to go upward at 45° 
except for very large horizontal gaps, where the concentration of 
fields at the side of the centre conductor causes them to level out. 


2 = 0.125 contour. 


Some scatter is visible on the & 

The attenuation is proportional to the product of £e and 
F (>) so this is plotted in Figure 3.9. There are some dramatic 
changes in direction of these contours, and the product is less 
than unity. (If F (Me were unity, it would imply that all the 
power is carried beside the centre conductor.) This plot, while it 
tells the designer how to make an applicator with a desired loss, 
does not indicate how to choose among the infinite number of 
possible combinations on one of the contours. 

In Figure 3.10 the contours of constant characteristic 
impedance (Z,) are mapped. Along one of these contours are various 
applicators which can be cascaded without reflections. These 
contours cut at a large angle across the contours of Figure 3.9, 


and thus permit a more intelligent choice of gap widths. In the 
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POWER CONCENTRATION FILLING FACTOR PRODUCT 
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- Using this map, the reflection coefficient may be calculated. 
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region to the right of the diagonal through the plot, the upper 
gap dominates the characteristic impedance and the designer can 
choose the horizontal gap for attenuation without affecting the 
characteristic impedance significantly. 

Since the applicator should be reflectionless, the 
interpolating program also prepared data for a and the F (Me 
product along the Zo contours. In the next two figures, it is 
plotted. The independent variable is taken to be the horizontal 
gap and the vertical gap may be determined from Figure 3.10. 

In Figure 3.11 the ordinate is ae and in Figure 3.12 it is the 
F(M)e* product. 

Again, ae can increase without limit while the Fee 

contours are limited to unity. In their middle regions the ze 


é curves behave as l/d. For 


curves behave as 1/d? while the Fé 
small d the departure of the Fee curves from 1/d behavior is due 
to the departure of the Ly contours from constant upper gap c. For 
large d, the fringing fields near the side of the centre conductor 
carry more power than the rest of the web region and the side wall 
of the outer conductor does not appreciably affect these fields. 
Thus there is a minimum attenuation which decreases as the upper gap 
is decreased, giving reduced Zo: In Chapter 2, it was shown that 
the rectangular hyperbola region gives this applicator its advantage. 
These graphs show that the low impedance contours give the designer 
the widest range of attenuation control. 

To use Figure 3.12 with constant upper gap c, a dotted line 
has been drawn across the chart showing the limit of the low 


reflection region. This line corresponds to the diagonal in 
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Horizontal gap, d 
Figure 3.12 Chart for Web Power Absorption Factor - no web gap 


- The dotted line is the boundary below which the upper gap may 
be made relatively constant. 
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Figure 3.10 and along it the characteristic impedance will be 
48% too low, giving -9 db reflected power. The chart errors are 
minimal except in extremely wide applicators, where the matrices 
were ill-conditioned. This error jecvioiiis along the 14 ohm 
contour. At the departures from hyperbolic behavior at low 
attenuation the errors are very small and the useable ranges are 
accurately given. 

To test the accuracy of the main program, a short 
section of rectangular coaxial line was constructed. Using an L-C 
meter and a time domain reflectometer, its characteristic 
impedance was found to be 46.6 ohms + 0.5%, when the side and 
upper gaps were equal. The main program predicted 46.5 ohms, 


(48) As the side gap was adjusted 


which agrees with Metcalf 
from 1 to 4 times the upper gap, 2 to 1 attenuation control was 
observed, which agrees very closely with Figure 3.12. 

The upper and lower gaps are normally used to control 
the structure, but if the web gap is large, the electric fields 
to the side are reduced. To check for this, the upper gap was 
fixed while a variety of web gaps were used in the main program. 
This results in applicators which do not follow a Z, contour 
over the entire adjustment range, but do so if this range is 
kept below the limit shown in Figure 3.12. In Figure 3.13 the 
upper gap was set to 0.2 of the centre conductor and the web gap 
was 0, 1/10, 1/4, 1/3, and 9/10 of the centre conductor. The 
ordinate is the F (Me product. Five to one attenuation control 
can be attained with very small web gaps but the presence of even 


moderately wide web gaps seriously degrades this control. The 
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Figure 3.13 Effect of Web Gap on Power Absorption 


-Upper Gap = 0.2m,about 25 ohms. Web gap = 0, 1/10, 1/4, 1/3, 
0.9m, centre conductor assumed equal to Im. 
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lower attenuation limit is also reduced by very wide gaps, a fact that 


can be used to further reduce the attenuation. 


3.4 Design Example 


To construct a simple applicator which will accept a wide 
variety of loads, several passes should be made over the web at 
increasing attenuation. This will permit uniform heating at first, 
with gradual increase of attenuation when the energy is partly 
absorbed. The efficiency will be kept high by the increased attenu- 
ation. 

If a 3/16" gap is used for web transport, Figure 3.13 
shows that the centre conductor should be at least 2" wide. From 
Figure 3.10, the 20 ohm structure has a normalized gap of about 
0.17 above a unit centre conductor, or 0.33" spacing to the upper 
ground plane. Checking Figure 3.12, little is gained from side 
spacings greater than 2 units or 4". By staying below the dashed 
line in this figure, we may avoid varying the upper spacing, which 
would otherwise be needed to maintain a constant Z. However, we 
lose a factor of 2 in attenuation control, and a pass will be 
included with a smaller side gap than the vertical gap. To gain 
a further factor of 2 in the first pass, we may add one section 
with a very large gap for the web. This also reduces the leakage. 

The resulting applicator is shown in Figure 3.14. The 
centre conductor is 2" square and, except for the last pass of the 
coaxial line, the upper gap is 1/3". On the first pass, the side 
gaps are 4" and the web gap is 1 1/2". On the second pass the web 


gap is reduced to 3/16", doubling the attenuation. On the third 
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pass the side gap is reduced to 1" and on the fourth and fifth pass 
to 1/3". On the fifth pass the upper gap is increased to 1"*. 

The applicator designed here is only 26" long and 4" high. Its 
width need only be about 4" wider than the web to be dried. 

The rest of the design consists of choosing the water 
load, matching transformer and the lips to prevent leakage. Air 
must be provided and it is probably advisable to run 1/2" poly- 
ethylene tubing inside the centre conductor rather than simply 
blowing in air. This tubing has the advantage of supplying air 
more evenly along the length of the web, and holes in the tubing 
will provide jets which impinge vertically on the web rather than 
flowing over it in a laminar pattern. This could also provide for 
air suspension of the web. 

In the low attenuation passes of this design, the main 
program: predicts 2000 volts RMS per cm at 10kW. Since the mesh 
used was as coarse as that shown in Figures 3.5 and 3.6, this 
figure is optimistic, and applies to a radius of curvature of about 
1/4" on the corners. Allowing for double the predicted field, 
j.e. 1/8" radius of curvature, the structure could probably carry 
2 1/2kW with 2000 V/cm RMS field strength. Circular 50 ohm 
1 5/8" coaxial line can also carry 2 1/2kW at 2000 V/cm RMS field 
strength. S-band waveguide can carry over 100 kW of power at the 


same RMS field strength and L-band waveguide can carry over 1 MW. 





* There are two reflections at the ends of the square pass 
which will cause up to -6 db of reflected power, depending on 
frequency. This has been tolerated to simplify the design. Also; 
in the next chapter it will be shown that using two passes at each 
attenuation sharply reduces inhomogeneity of heating, and more than 
One mode can exist at 2450 MHz. 
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CHAPTER FOUR 
INHOMOGENEITIES AND ATTENUATION IN TRAVELLING WAVE APPLICATORS 

In Chapter 2 it was shown that attenuation could be control- 
led in TEM structures and in Chapter 3 detailed design data were pro- 
duced for this purpose. In this chapter information will be given 
to select the optimum attenuation. 

An electromagnetic field can give uniform heating of homog- 
eneous webs if there is no reflected wave to cause interference patterns 
and if the energy in the incident wave does not become absorbed too 
near the edge of the web. In addition, if a web has wet spots they 
absorb more power than the rest of the web. This fact permits levelling 


of the moisture content of sheets such as mapertwaen 


Again the field 
must be kept uniform. When large reflection coefficients occur in a 
travelling wave applicator, or when the web attenuates the wave too 
rapidly, the result will be poor product quality and possibly fires. 
A more accurate technique is needed to predict inhomogeneous heating 
in wet webs. 

Heenan! 392 treated the situation where the evaporating web 
moves in the same direction (concurrent) or in the opposite direction 
(countercurrent) that the microwave power travels. He did not go 
into detail about temperature, vapor pressure and air flow; or about 
reflections. To include the latter causes of inhomogeneities, it 
will be desirable that the web be stationary. This is batch heating, 
which is not used in practice, but which gives a worst-case inhomog- 
eneity prediction since the motion of the web will even out the heat- 


ing. If the predicted inhomogeneity is low, Heenan's theory can be 


adapted to crosscurent motion by considering slow waves. 





-fowtnos ad biuoo nofisyunotis deste rworke 2W st s $ nad a 
Sate’ ‘@ i? 


-o1g 9t9w $i8b nptesb belrsteb € yedgsdd nt bas sorus.unte er : oe ae 
navip sd Tliw nottemotat yetqedo etdd ol .seoqiwg tag not bsaub 
-Nofteunstts mumitgo at tosfse of 

-pomod to paftsed mrottnw avip aso biott sitsnpamoytosis nA 
envettsq sonevetistnt seuso ot ovaw befosltoey on af stort ld hanes cyoans 
oot badvoeds smoosd ton esob svew dnsbtont odd nf ypyvens on? 7 oe 
yond 2toge tow 2s dow 6 ti .mofttbbs al .dew edt to spbe ont bien 
pniffevsl 2dtmisq fost efdT .dew ond to test ants neay "reMaR s710Mi dno2ds 
bfoit sit mispA EM sang 26 flove efsend2 to Ingtnoo syutaton ony 0 
6 Nf Wwo00 esnatortts09 notsasltsy sprsf nenW .myothau tqod od teum 
oot svsw sit estaunotis dow ot nadw vo engines eee Svew enti fevers 
-esttt yidiezoq bré ystlsup toubovq 100q sd TT tw + huean ong — 


pnitssd evosnspomorat totbearq of babosn 2t supindoas sil ci gvom A 
: 7 Wee & 


daw Ont sIeyvegsve eds 9ovorw NOFI HOT te ont batsend (Ones 
nofsisextb st 20qqo eit mt vo (3ns9w2N09) notsoarttb ame ont nt avon 
op don brb oH .2Tayens ‘3W0q eveworatm ond Sand (nerorstaues) 


_ tuods yo iwolt vis -bis hate von aim mys act 





63 
4.] Drying of Wet Webs 


In conventional hot air drying there are two stages of 
moisture removal, called constant rate drying and falling rate drying. 
The moisture is removed first from the surface of the web, and capil- 
lary action keeps the surface supplied with moisture even though the 
density of the web is falling. During this constant rate phase, heat 
is supplied directly to the wet surface and it is usually carried by 
the same air which carries away the water vapor. Later, the evapor- 
ation is carried out inside the material and heat is supplied more 
slowly to the water. This is because the heat must penetrate a dry 
thermally insulating layer. This is called the falling rate stage. 

Recause energy is absorbed from the electromagnetic fields, 
mainly by the water itself, migration of the water in the material 
has little effect on the power absorbed. If the power is travelling 
forward only in a travelling wave dryer, the evaporation rate begins 
falling as soon as drying begins . However, the rate of diffusion of 
water vapor inside the web will limit the transport of vapor, result- 
ing in a higher temperature. Thus a wet surface is advantageous, 

but not as much so as in hot air drying. 

The moisture content above which free liquid water is mobile 

in a web is called the fibre saturation point 5 Above this point 
ef. of the water component of a web is approximately the same as for 
free water, while below it the large percentage of bound water causes 


increased Ep isto Bia Phe For wet webs, high temperatures reduce 


a 


*After transient conditions disappear: see later. 


sFor Douglas fir, below about 50% dry weight or 0.2 gm water per ml. 
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a while for drier webs high temperatures can increase it (Ol), This 
means that there is the danger of thermal runaway and fires if dry 
spots develop, and the applicator should be designed to prevent in- 
homogeneities. 

While the constant rate drying stage is not present in trav- 
elling wave microwave driers, the same range of moisture content per- 
mits the use of two simplifying assumptions; namely that the required 
vapor pressure (and therefore temperature) is determined by surface 
conditions and that the lossiness of the web is determined by free 
water . Below fibre saturation these assumptions break down, and the 
required vapor pressure rises significantly at about the same time as 
Ey does. The analysis of this chapter will predict inhomogeneities 
in this range, thus determining whether there are spots which have 
dangerously low moisture content while other parts of the web need 
additional drying. 

Another analysis of inhomogeneity® assumes that the web is 
uniform, that there are no changes in temperature, that no reflections 
occur and particularly that no evaporation takes place. This analysis 
will be derived in Section 4.3 by using simplifications of the more 


complete treatment. 





*If a multimode cavity is used, the electric field can be increased 


(31) 


sufficiently to give constant rate drying This occurs automat- 


ically, even below the fibre saturation point, if the efficiency is 
kept high and the magnetron power output is constant. 
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It would be possible to relate vapor preceunen power 


absorption per unit length, dielectric properties of water'20) 9497) 


(68) | and other parameters of a specific system to predict its com- 
plete behavior while heating a wet, moving web. However, in the fol- 
lowing section it is assumed that the web is batch heated, and to 
further reduce the complexity, the remaining variables can be normal- 
ized to give only two independent variables; the initial equilib- 
rium temperature and a product of filling factor, moisture content 
and field concentration factors. The result is a moisture contour 
which is initially level, becomes uneven as the edge of the web near- 


est to the magnetron is dried, and is eventually levelled out when 


the energy passes further over dried parts of the web. 


4.2 Energy Deposition and Drying for Very Wet Webs 


(a) Power Applied 


The power dissipated in the web per unit length is 


*n dS (4.3) 


and 
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. . . 2 
which gives the field concentration due to reflections. & is 


related to fields and power as in Chapter 3. The wave impedance is 


E uo/e, 
n= f/He = ———_— ory /ev 1 -(w./w)? 4.5) 
ei | 


for TE and TM modes respectively, assuming that very little power 

is carried in the web. The complex dielectric constant of the water 

in the web is 
e = e fey, - jwe,,) 

and p is the volume fraction of the web which is filled with water. 

The a of the web as a whole is pey. We is the radian cutoff frequen- 


cy and w is that of the power source. 


(b) Drying Rate 
When the web surface is kept soaked by capillary action 


the mass lost through drying is given by the empirical relation 


Q 


e = Ke (R)"(s)" (Pa/pe,) py B (4.6) 


Q 
c+ 


*In this analysis it is assumed that the power absorption is caused 
only by free water in the web. Thus Fo is the cross-sectional area 
occupied by water with dielectric constant oN ~ Jegey: Alternat- 
ively F is the filling factor of the web and its dielectric constant 
is e,(e,- jes). The latter interpretation may permit easier under- 
standing of the problem at low moisture content, where ec is higher 


than predicted by this model, since ey is not constant ‘6 ; 


(Ba) 
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where m is the mass of water per unit length, K, is the mass trans- 
fer coefficient, R is the modified Reynolds number, S is the Schmitt 
number, B is the boundary exposed to air, u and v are experimentally 
determined constants, Po is the total pressure of air in the (dry) 
air stream, Py is the partial pressure of water at the web boundary 
and Pen is the logarithmic mean partial pressure of dry air in the 


boundary lamar’. The variable Pen is given by 


= (Peg = Ppt) (4.7) 
an (De./Pey) 


Pm 


where Pre is the partial pressure of dry air at the surface. Pee is 


the partial pressure of dry air at the upper edge of the boundary 


(1). 


layer. Simplifying and combining drying parameters, 
ap gir B Pa Py (4.8) 
ot ASG, Pier 

where Kary = K.(R)"(s)", which varies with air duct design and flow 


rate, Ww is the density of water (kg/m), AS is the cross-sectional 


area and B/AS is the surface to volume ratio of the web. 


fe) Temperature Changes 


The rate of temperature rise is given by 


Be cats, eee Ns 9p 
ot pASG C,, at pASG.C,, (Pureb a wwe at) (4.9) 


where 6 is the average temperature (°C), Cw is the specific heat of 
water (Joules/°C kg), Q is the heat in the web (Joules/m), r,, is the 


latent heat of evaporation (Joules/kg), where MKS units have been 


v9 
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used for consistency with the electromagnetic units. It is assumed 
that no heat is lost to radiation and that none is supplied by 


the air stream. k is nearly independent of Lonivarature lis 


dry 

Let the equilibrium temperature be 6, If 6 is below 8° 
de/dt is low and p changes slowly while 6 rises. If Poeis is constant 
and Ae = (6 - 6.) << 6,5 a change in dp/at will occur, i.e., 


phoc, = 2 f° {{28} - {22} ,} dt. (4.10) 


Differentiating with respect to time, 


a0 = Ww 


90 _ (90) 
at Cc ( ( 


4.11) 
ate) 


Thus the time constant of temperature equilibration is 539.5 times 

less than that of drying of the web. It can therefore be assumed 

that the web is very close to thermal equilibrium while it dries. 
The vapor pressure of water and Ey are known functions of 


(26) (57) ,(68),(71) Since the surface is assumed to be 


temperature 
wet, Pv is the vapor pressure, and from a curve fitting program, 


6 


p, 2 0805x1078 95 + 2.097x10°© 6% + 1.942 10°* 63 


+ 0.011084 62 + 0.32593 6 + 4.6249. (4.12) 


For the dielectric losses, at the 2450 MHz ISM frequency, 


-] 


Ew + 3.6949x10 ° 64% + 0.0001084 63 + 


0.01222 9+ = 0.70285 8 + 21.331 . (4.13) 


The dimensions are, of course kept correct by the numerical 


constants. 
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(d) Normalized Drying Rate 
Substituting 96/at ~ O and using Equations 4.2, 4.8 and 4.9, 


Wy 22 é p 
E Pe wyrer n Pe F aa Kary B ( a/Pe) pe (4.14) 
If the web has a thickness t'and two sides are exposed to air, 
BAAS = (2/%: (4.15) 


and we may define 


b, (0)? Pye) (4.16) 


where 


Uy = 1 Torr (Az) 


which keeps the dimensions correct. Then solving Equation 4.14, 


2 Wowk AS Me 


Sykes p dry (4.18) 
e t Egnw Pe WoPE 


The variable xX, 1S dimensionless and depends only on 6: Solving for 
X,at Pe. = (760 - Py) and a 760 Torr, and changing dependent 


to independent variable with the curve fitting computer program, 


Pos -0.01800261 x* -0.08522502 x3 +.4938560 x? 


-10.37308 x +15.19636 (4.19) 
where 
x = on(x/pE*y2(Pe/P)) (4.20) 
and where 
X(8,)- (A520) 


This gives two new variables, P_, a reference power, and 84° the 


6) 


reference equilibrium temperature at this power for p = y* = 1. 
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The other independent variable is Kary? which can be adjusted by 
changing the air flow. While it does not appear in the following 


analysis, it will be assumed that the air flow gives an equilibrium 


temperature 8. at the forward power Pos From Eq, 4.14, at equilibrium, 


Of = 6 hem ue E° 7 Pe PR 
at Sa ky ee (4.22) 
wowss 
Let 
ee 
= r (4.23) 
Tg 
where 
T A,.g, AS 
O " 
Ww 'T g nPwel(o )E2F (4.24) 
Then 
PoP ee?) 
1\/a09 f “wie 
(—) (=) = - y? (4.25) 
OT 0 ey (8,) 
eo Pe : 
i. represents the boil-dry time constant at ier) = y2 = ple in 
O fe) 


Returning to the temperature rise equation, 4.9, we can 


normalize it also: 


ge Pa bd Ew) ap (4.26) 
oT P eh oT 
(6) wW Ey 0 
where " 
(Y) ee (4.27) 
T = — w i ° 
6 Cy Ew ae 


This can be used for the transient problem, if needed, and again 


shows that change in 6 is much more rapid than change ino. 


* The substrate is assumed to be lossless. 
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(e) Propagation Constant 
Altman (2) gives a perturbation solution for the propagation 
constant. Assume that all the web in a cross-section is at the same 
temperature and moisture content (batch processing). Then for a 
single forward or reverse component E. or Eres the components of 


the propagation constant are 


Eg ley -1 )we jee, SE) ds 


> So (4.28) 
2 [J (ExxH)*n ds 


and ff er 
€ € "wo (E*-E) dS 
_ SOS Bee (4.29) 
2 Jfg(E*xH)-n ds 


where n is a unit normal vector. EG is the unperturbed electric 


field while E is the perturbed field. For lightly loaded LSE modes 


EG is approximately equal to E. Also, E = nHxn and a wl! = 


tee i / vi - (w/w)? . The propagation constant is then 


W a ' F 
ce" nune2py? + rege? e,le)-1) nwxt*pv?) 


P (4.30) 


= + 7 = 
i esa bef: Eo 


where Vp is the phase velocity in the waveguide. 


[) Capacitance and Characteristic Impedance 


The capacitance of a transmission line can be evaluated 


by using energy storage, or (30) 
Si as 
meehis_t- EF [J (vo-vo) dS (4731) 
Vv Vv 


where v is the total voltage between conductors, No is the stored 
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electric energy per unit length and C is the capacitance per unit 
length. For a perturbed TEM line 

GC lly BxD! Jf (EX +E) ds 

Ci = (ey, - ] Diseeras aa A Ay 
—0 k. 
Jf, (es E) dS 





HE (ehd-T)e? Fo (4.32) 


This can be used to calculate the characteristic impedance change, 


Z, =vL/C, and.Z/ = vb (GE, Ac), (4.33) 


In complex form the normalized characteristic impedance is 


Z 1 
£=T 7 VV + Fetes - 1 - je") (4.34) 


(g) Distributed Reflections 


(14) that a TEM line or waveguide with 


It can be shown 
a characteristic impedance discontinuity and a load mismatch (Fig. 
4.1) gives a voltage reflection coefficient 
r, +r3exp(-2y2) 
r= —Trexpt -ayn)* T; + T3 exp(-2yz) (4.35) 


For a transmission line with gradually changing ZW? we may assume 


step increments at the ends of infinitesimal lengths dz, or 


] 
a d gn Zo (4.36) 
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or at the input of the line or waveguide, 
3 


F 5 Tp exp(-2yz) CS an a dz (4.37) 


ae! 
" 


Using Z, Eq. 4.37 may be approximated by two finite difference 


equations in Ep and E 


f° 





Figure 4.1 Double Discontinuity 


(h) Finite Difference Equations for Reflections 


Assume that the applicator is matched except for the effect 
of the web, ie once the power passes the load it is absorbed. The 
web may be divided into N segments, h units in length, in the propag- 
ation direction, and numbering will begin at the end opposite to 


the power source. 











% 
re el emer iies Ke8 
EY Ue a : Er, ey al 

ie . 
N a m-t ie t 


Figure 4.2 Distributed Reflections 
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At each segment, assume the reflection to occur at the load 
end and measure E. and E. at the beginning. Using Eq. 4.34, 
My (4.38) 


en eae: 
Starting at n = 1 we assign E. some value temporarily and calc- 
] 
ulate all the E. and E. by recursion. (Finally, all values can be 
n n . 
reset to correspond to the actual input power.) Then 


ER = ERT exp(-2y7h) (4.39) 


For any element n, from n = 1 to n = N-2. 


E 
phetros 1 ey AMA rate suph, (4.40) 
n+] ce a ose ee beeen 2 
Clit ors) 
n 
and 
E TNE rT. ,,exp(-2y.h) + E. (1 - r_,,)exp(-y_h) 
aE tom fat] n+] n es n+2 n 


(4.41) 


For the last element, Equations 4.40 and 4.41 give the correct res- 
ults if n+] = N and Tus 0. N, of course gives the web width,z_, 


j.e. N= z/hs to the nearest integer. Then 


2s 
Yh ~|—2—_®) -|_2__4 (4.42) 


(i) Finite Difference Time Response of the System 


Eq. 4.25 can now be converted to a finite difference equation 


by setting E, Ee = ] for normalization of power. A small inter- 
N N 
val of normalized time t can be chosen. Then 
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w 
- 2 k,n 2 
(KT 7 Pn =P kn en (60) ve kan (4.43) 
At 
where the subscript k refers to the time index and n refers to the 


Zz index. 


At k = 0, 84 and some desired Py are assigned as initial 


conditions. The web-applicator-magnetron system design is described 


by the constants 3) hh 0 > W.> Wy, P_ and the air flow. The 


0° 7m? %¢ 0 

web reaches an initial equilibrium temperature of 8, at n= N, or 
where the power is first incident. The air flow and EO are 
assumed to be correctly related to give 85° Then the edge nearest 


sl sh 


to the magnetron will dry at a rate such that (0, nl a 
N = SN 


where the subscripts o and N indicate initial conditions at the edge. 
The air flow is assumed to be uniform and initially dry, 

the web to be thin enough that it is at nearly the same temperature 

throughout, which gives the vapor pressure of the surface water, and 


k is assumed to nearly constant. Only the range above fibre 


dry 
saturation is treated and diffusion of liquid water from location 
to location is ignored. Perturbation analysis of waveguide behavior 


has been used with LSE modes. 


4.3 Exponential Heating Model and Meanderline Applicators 


There are important cases where evaporation is minimal or 


irrelevant. Then Eq. 4.9 reduces to 


78.3 ] 
all ASG uC 4 ree (Leff) 


where 94 and Cy are the mass density and specific heat of the web. 


(eh.5) 
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Expanding P 


—web’ 
10). .0 6 (8) we Fy°P, . ’ 
e) ASGyCq (4.45) 


In this section a (9) is that of the web as a whole, and must be 
obtained from sources such as Tingal ©!) Integrating Eq. 4.45, 


ue Fue t " 
6(t) = ASGuCy Je. Pe Ep (6(t))dt 7 (4.46) 


Should “ be nearly constant, 


Ege (EFF) 
6(t) = The Pe (t-t,) (4.47) 


where mq is the mass in the applicator per meter. 


The wave is attenuated in propagating along the applicator 


2 


according to Eq.4.32. Assume ~» = 1; i.e. there are no reflections. 


Then, as in Chapter 1, 


in =) Pe exp ( -2az) 
iW fla 
ca exp (-e,e,nw (Fe-p))z (4.48) 


where i is the power produced by the magnetron. The temperature 


profile as a function of time and distance is 


PG(t) 
e(t,z) = Maca {t-t,)} exp {G(t)nz} (4.49) 
where G(t) = eae phe) Nee is (4.50) 


The variable G is the conductance per unit length and is the conduc- 
tance of the web multiplied by FES Gn is the ratio between power 
absorbed per unit length and the transmitted power for the system 


under consideration. 
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When many passes are made, for a web Zn meters wide and 
for a residence time t, at location z on pass k, 


PG, (t-t,) 


At) or rene exp|-Ginz] + exp[-@z,-z)G,n] 


+ exp [-(2z,+z)G4n] + exp[-(4z, -z)Gyn] re ek (4.51) 


The temperature inhomogeneity is less than that caused by a single 
pass, which is also true for moisture content inhomogeneities. If 
the series Gy (G).G, etc.) is gradually increasing, the range of 
loads which may be heated uniformly is increased. 
In a Taylor series expansion for 6 the first order term 
may be cancelled. When all the G's are equal this can be accomp1ished 


(73) 


by using an even number of passes , and the decrease in inhomo- 
geneity is dramatic if the total inhomogeneity per pass is small. 

For instance, 20% inhomogeneity can be reduced to 5%, 30% to 10%, 
and 40% to 26%. For such an applicator, the designer should check 
for excessive reflection by using Eq. 4.34, check for excessive 
inhomogeneity by using Eq. 4.49 or the drying profile of Section 4.2, 
and if these do not exceed about 30%, add a final pass whose attenu- 
ation is sufficient to cancel the linear term of a Taylor expansion 
of Eq. 4.51. Alternatively he may use an even number of passes at 
each attenuation rate. If the inhomogeneity is excessive for any 


pass, it may be reduced by reducing the product of F and = for the 


particular applicator chosen. 
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Figure 4.3 shows the percentage temperature inhomogeneity 
(lower three curves) and the percentage efficiency (upper three 


curves) for three applicators. The inhomogeniety is given by 


(4.52) 


= 100 Beat - 6. )/6 


min 
and is plotted on a log saale The abscissa is Gnz, also on a 

log scale. Numbering from the top, curve number 3 gives both 
efficiency and inhomogeneity for applicator Al» a single pass of 
waveguide, TEM or any other mode. The user is given the choice of 

bad inhomogeneity or bad efficiency which is not surprising. 

Curve 1 is the efficiency and curve 4 the inhomogeneity for applicator 
Ao consisting of five passes of the same waveguide as in curve 3. 

This applicator gives improved uniformity since it returns the wave 

to the far side of the web, and improved efficiency due to the 
increased number of passes over the web. Curve 2 gives the effi- 
ciency and curve 5 the inhomogeneity for the graded attenuation 
applicator, Az, designed in Chapter 3. The abscissa is Genz s where 
G,; Go» Gs, Gy» and G, are increasing in the order of 1/16, 1/8, 1/4, 
1/2, and 1. There is some loss in efficiency because of the low 
attenuation on “ae passes, but there is a more than proportional 
improvement in homogeneity. The final pass to cancel the first- 

order Taylor series term has been omitted. It would produce better 
cancellation in the graded applicator A3 than in the constant 


attenuation applicator Ao. 


4.4. Experimental Results 


A length of WR284 waveguide (1.34" by 2.84" inside dimensions) 


was constructed with z wide slots in the broad walls and tapering at 


8t 
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the ends, with a useable length of 45 cm. A water load, matched to 
provide less than -20dB return loss, and a Philips Dx260 1 kW 2450 MHz 
magnetron were connected to either end of this TE, applicator. An 
air duct was provided with an opening into one slot x high by 18" 
long. Five vanes were mounted internally to redirect an air stream 

1 la 


> by 15 end evenly over the web. This was checked 


with a simple flowmeter. The web was placed on a 0.01" mylar strip 


incident on the 1 


to expose only one side to the air stream. A Fortran IV program* 
was used to predict the attenuation and temperature profiles given 
the initial temperature, which was measured after 5 seconds (rt=0.2). 
Three strips of blotting paper were soaked to 185% water 
content (dry basis) and placed on the mylar strip. One kW of power 
was applied giving an initial temperature of 45°C where the power 
was first incident on the load. The experimental reflection coeffi- 
cient, measured at low power, was 3.2 +0.2 dB, or 45% reflected power. 
The program predicted 5 dB, or 32%. This is probably due to higher 
order modes being generated in the vicinity of the interface, since 
the load was a large perturbation. The predicted es was 17.8 sec at 
1 kW and the observed value was between 185 to 215 seconds. In 
Table 4.1 the weights of samples from the partially dry web are 
compared to the theoretical values. These results check very well. 
Due to the drying rate and the available balance used, the webs were 
heavy and measurement times in the vicinity of t = 0.75 were used. 
The moisture content is expressed as the percentage of the moisture 
content measured in an identically prepared control sample. This is 


of the same form as the output of the program, which assumed that 


the web started at p = 100%. 
* See the following section. 
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Table 4.1 Experimental Check on Attenuation Inhomogenéity 


Moisture Content 


. * 
Time Temperature Location Experimental Theoretical | + 


15sec _35,0° 2 cm 42+5% 39.3% 0:8 
9 cm 86% 91.9% 
5 18 cm 97% 100% 
15 sec 35.6 1 cm 45% 44 8% 0.7 
250 5 cm 73% T6215 
9.99 9 cm 84% 93.2% 
5° 14 cm 97% 99. 9% 
50 30 cm 94% 100% 
50 39 cm 104% 100% 
‘ Measured temperature was 35° at/r = 0. 


4.5 Computer Results for Very Wet Webs 


Fortran IV programs were written to solve the system of 


difference equations, developed in Section 4.2, for TEM and WR284 


applicators. The input variables were ec F, the number 


Ww 
of wavelengths of web, the time and distance intervals, and printing 


(05) Xo» 
instructions. The output for each time interval was maximum 
reflection coefficient and, at several locations, the temperature, 
water remaining, and deviation from an exponential which was fitted 
to the curve of moisture content. 

Figure 4.4 shows the moisture content curves for one wave- 
length of wet web at several instants in normalized time. Feo was 
0.03416, giving Fe*c, = 0.1 at 80°C. +, starting at the top, is 
Osea Op5 Ss POG ee Lehoe Veoe 1.7) ANG oO 1Or che 90) Tecat Ves 
representing g0°c. The x's are for t = 1.0 at 35°C. The frequency 
was 2450 MHz and the applicator was a TEM structure. Waveguides 


give more attenuation in the TE modes. 
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Moisture Content at Time 0.25, 0.5, 0.75, 1.0, etc. 


The moisture content starts at 1.00 throughout the web and 

evaporates. Plots are shown at the instants indicated, and 
the hot spots are lower in moisture content. The wet spot 

at 1 meter has been shielded by reflections and appears in 

all plots. 









a Saslninlid- ends dr abitin Onetane ewwsetom 


yp cendgommacbed pwolz 218 2of4 .2e36-10qsvs 





-dnsdnos swuserom at teol 94s re * a 7 


83 


One rather bad wet spot remains at the far edge. Since a wet 
spot is a region of low electric field, it is shielded by reflections. 
The fact that the guide wavelength increases as the water evaporates 
prevents standing waves from increasing the shielding effect as evapor- 
ation progresses. The maxima move, drying the wet spots previously 
left by minima. This effect would not be present in moving webs. 

Other results for webs wider than 1 wavelength show that at 
2450 MHz, inhomogeneity is introduced mainly by attenuation. This 
shows that at least in the 2450 MHz ISM (Industrial, Scientific and 
Medical) band, problems with standing waves are likely to originate 
in the design of the unloaded applicator, not in the reflections 
from web discontinuities, as had been suspected. The reflection for 
the edge of the web in Figure 4.4 is only -11.6 dB, which is not 
difficult to create by designing a bend in a waveguide. -11.6 dB 
is not a severe reflection, but its effect is quite substantial due 
to the fact that the VSWR is twice the magnitude of the reflected 
electric field and the power absorbed is proportional to the square 


-of the sum of the reflected and incident fields(72), 


The upper gap 
controls the characteristic impedance in the TEM applicator, which 
should help reduce reflections. 

The least-squares fit of an exponential through the moisture 
content curve of the web showed that the heating can deviate marked] y 
from an exponential due to high temperature depression of cy: As a 
result, the drying is more uniform than the heating predicted (Section 
4.3). If the air flow is increased to keep 6, low, levelling is 


improved, but an uneven electric field will then create more 


inhomogeneity. 





















8 | as i. 
tow 6 sont2 .opba vet ond $6 antemay soqe tow bed wane ey, “vi an 
_enotigettss vd bebfonda at tf .bISt? ofttools wolf Yo notpsy 6 ef doge 

2etsiogsve vetaw ont 26 292n970Nt diensfoven ebtup edd ents tos? sdT 
-oghve 24 Soette pribiatie odd pntesstent movt eovaw onthaaty 2tneverq 
yfevotve iq etoge tow sit pratytb .ovom smixsm sAT _eszeatporg Notts 
.2dew pnivom nf tases%q sd ton bluow Sostte 2tdT .emintm yd Stef 
$6 tend word? dipnefsvew [ nett vobtw edaw wot atfuesy r9eds0 
2inT = .nofisunetts yd ylatem beoubovint 2t ytrenapomorat eSHM O@bS 
bis oFtrineto2 .[6tvs2ubal) M2T' sHM O@bS edt at tesol 36 Jans zwone 
atsnipiyo di yfextl os 2evew ontbnste dttw emefdorg . basa (Teotbom 
anotiveftet odd at Jon ,torsdifqqs bebsofau art. to.aptesb ont at | 
0? Noffosftey eiT .betosgaye need bed 26 sestdtunttnooeth dew mort 
ton et dotiw .8b a.ff- yfno at &.8 swuptt nt dow oft to spbs orig 
fb 3.1f- .sbhupsvew o mf bred 6 pntnptesb yd sens ot dTyottttb 
uh Tsidnstedue edtyp ef tse%te ett sud .not}oefter sisvse 6 ton ef — a) 
batosltey oft to shuttapsm edt eotwt et AWeV ols Sard 3267 ont of 
sTaupe sit oF Fenotsroqorg 2t badtoeds T9H0q edt bis biatt diisoats 
qsp veqqu odT (EY) sp rart tnebront’ bas bstoslts1 sdt to mue oft to: . 
dotiw Nossai Faqs M3T oft nt sonsbeqnt oftetratosyendo ont elovtees | 
_ - s@notsaafte soubs1 ofed Ohaate z 
eibiehotbiteardts Svadeans n6-10 3tt 2sysupe-desel oT oa! 
ab 69 Git ded as “Ferd Bavodd daw aff 'to avo 9 
i. hes seh. ow suserogmat Ae td easter 


a tn or 


arene | ~~ @ Pease. 


e>4 













: 7 Sheena > 2 : ; 7 
; st oe ae) a 2 
¥ : 7 - bs i 





84 


In Figure 4.5 some results are shown for inhomogeneity vs. 


normalized time. The curves are: curve 1, Fees, (e)= 0.03416 or 


Fe aor (0) = 0.1; curve 4, with x'es Fee = 1.01113 or Feée, (6 )= 


6) 
0.03333, where the waveguide version of the program was used; 


2 


curve 2, also Fe =0.01113 etc. but for TEM wavequide; curve 3, 


3, Fee = 0.003416 or Fee (0,) = 0.01. In curves 1, 3 and 4 
a TEM version of the program was used and in curve 2, the waveguide 
version for WR284 guide at 2450 MHz was used. The initial temper- 
ature was 80°C, and 3.33 free-space wavelengths of web were assumed. 
The ordinate is the difference between the wettest and driest spot 
in the web, referred to the original moisture content, or 

fd (max f 


ty ts 
—— xX 100% (4.53) 


where p was initially 1.0. 

To use this chart for conditions not exactly as programmed, 
we note that in all but very short webs it was found that the slow- 
wave effect smoothed out the wet spots. When wider webs, in this 
case 10 free space wavelengths,were assumed the inhomogeneity was 
proportional to the web width. When different temperatures were 


assumed, the inhomogeneity was proportional to e, (6 ). Due to the 
higher n of the TEy waveguide, attenuation is approximately doubled 
for this particular guide at this frequency. Inhomogeneity is 
proportional to n/ng . Then for curves 1, 2 and 3, the factor Fee 2, 
is 0.333, 0.1, and 0.0333 respectively. The wet spot is not the same 
magnitude because the reflection coefficient is lower for wider 
lighter webs and some deviation will occur from scaling to narrow webs. 


We may even change the frequency, remembering that reflections may be 
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much worse in the 915 MHz ISM band, where c,, is lower. The TEM 
plots represent G = 2.1, 0.70, and 0.21 respectively. For curve 

3 of Figure 4.3 the respective temperature inhomogeneities for the 
other model would be 85%, 40% and 17%. Thus there is considerable 
gain in homogeneity from evaporation and the depression of e,, at 
high temperature. However, the latter effect tends to work against 
levelling of moisture. 

Figure 4.5 differs from Figure 4.3 in that the web will 
eventually dry and the power is assumed to go to the water load 
when it does. The unused power goes on to dry the remaining wet 
areas in this batch processing scheme. Batch processing is not 
often used in practice, but this assumption gives a conservative 
design plot. The average inhomogeneity in one pass will be 
slightly better for a moving web, since it will be partly wet and 
partly dried, i.e. different parts of the web represent different 
instants in normalized time. Therefore a bound for I in a moving 
web may be obtained by finding the appropriate curve in Figure 5.4 
and adding a term te The latter term is needed because in moving 
webs the standing wave nulls remain stationary after equilibrium 


is established. The extra term is given by 


Ripe a | (4.54) 


where T is calculated from p = 1 and p = O in Eq. 4.38. 
If the air flow is increased, the initial equilibrium 
temperature is reduced. This means that the effect of high temper- 


ature depression of Ey is less, and I is increased. The effect of 


mass transfer still exists, though, and a graph such as Figure 4.5 
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for the new temperature would be needed. An even more conservative 
bound could be obtained by using the L. in Figure 4.3, which does 
not allow for evaporation. 

In moving web systems one cannot always rely on eventual 
drying of heavy webs. Should the material be susceptible to thermal 
runaway, it is very undesirable that one part of the web be much 
drier than the rest. If the graduated multiple pass applicator is 
used, the web is forced to accept power by the very high aA in the 
last pass. If the power is too high and the material is susceptible 
to thermal runaway, fires may result. The web must not stop and it 
may be disirable to use limited attenuation. Increasing the air 
flow rate gives lower temperatures and may keep the web in a range 
where there is no danger of thermal runaway. This does not effect 
the drying time, except for an increase c, for the free water in the 
web. It is desirable that re be kept low, however, so that these 
effects do not occur. 

It is now possible to inte the evaporation analysis to 


(35) Several passes of waveguide form a slow 


Heenan's treatment 
wave structure. Let the distance measured along the wave path be 

z and that along the slow wave path be z'. The slow wave travels 

at a velocity which is reduced by Z'/Zys where Z'p and zy are the 
periods of the path of the fast and slow wave respectively, in 

other words the total length of the waveguide per pass and the spacing 
between passes. The wave travels in the +z and +z' directions. 


These coordinates are at right angles to each other. Heenan refers 


to any web motion which is perpendicular to z as crosscurrent flow. 
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For the slow wave, the material ard energy flows are concurrent if 
the web material moves in the +z' direction and they are countercurrent 
if the web moves in the -z' direction. 

Assume that the predicted inhomogeneity for batch heating 
is less than 30% (or some other accuracy criterion). Then an even 
number of identical passes will give a uniform evaporation profile in 
z. If the web is moving this profile will vary across the waveguide 
in the z' direction. Heenan's treatments of concurrent and counter- 
current flow are then approximately valid if z' is substituted for 
his z and all the attenuation factors are multiplied by 2/2). The 
applied power is a step function in z' and the moisture content is 
constant between passes and decreases across the waveguide portions 
of z'. The equations would predict an approximation to this profile, 
a continuous function equivalent to distributing the power uniformly 
over the entire web. 

In summary, the evaporation analysis has resulted in inhomo- 
geneity plots which differ by about a factor of 2 from the simple 
exponential. The observed hot spot formation at 2450 MHz has been 
traced to the applicator design, not loading, but at 915 MHz c, is 
significantly lower, and web reflections may be a much greater 
problem. The analysis can show whether this is so for any loading 
condition. It is not a complete description of drying below the 
fibre saturation point, where c, of water is not constant and the 
surface is not kept wet. It was shown that Figure 4.3 is a worst-case 
plot except in the range where de, /90 is positive, and that the 
inhomogeneity of drying can be conservatively predicted by calculating 


that in the initial heating profile. 
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CHAPTER FIVE 
CONTROL OF REFLECTIONS FROM RESONANT APPLICATORS 


The problems associated with resonant applicators will be 
considered in this and the next chapter. In travelling wave appli- 
cators, the techniques which were used for reducing heating nonuni- 
formity reduced reflections as well. In resonant applicators, 
however, the energy is reflected many times inside the cavity until 
it is dissipated. Conventional coupling techniques often permit 
some of the internally reflected waves to be coupled back out of 
the applicator. The resulting high VSWR shortens the life of the 
magnetron and is the most troublesome aspect of resonant applicator 
design. A coupling technique to reduce this difficulty will be 


treated, both theoretically and experimentally, in this chapter. 


bal Microwave Heating in Resonant Cavities 


When a load is relatively predictable, specialized appli- 
cators may be designed for it. However, experience has shown that 
for loads of various sizes and shapes multimode cavities often give 


(74), 


the most satisfactory performance In a well designed multimode 


cavity many resonant frequencies are available in the operating 


38) The electric field nulls of these modes do not in 


ISM band! 
general coincide, and heating will occur at almost any location if 
all the modes can be excited. Further, if a dielectric is placed 
in such a cavity, alternative modes are pulled downward into the 
frequency band of the magnetron, and the VSWR is better than for a 


single mode structure, which would simply be detuned. 
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The most common use of multimode cavities is in microwave 
ovens. Since these appliances are used by the general public, 
the system must be insensitive to moderate Soude W1>). Figure 5.1 
shows a typical loop coupled microwave oven cavity, and alternate 
coupling structures are shown in Figure 5.2. The filling factor is 
usually less than 10%. 

The loops and probes are usually less than 4/2 in size 
and do not extend over more than one field maximum. Therefore they 
are termed smal] element couplers here. The long strap coupler 
acts as a stripline transmission line with distributed electric 
and magnetic coupling. Because of its size it couples strongly to 


(75) 


more modes than the small element couplers and gives reduced 
VSWR. However, experiments show that the reflections are not 
reduced as much as might be expected and that stub matching in the 
transmission line is not uniformly effective over the entire band. 
For all the couplers shown reactive components of the input impedance 
vary widely with frequency; each change in load or operating frequency 
necessitates a change in the stub. This type of behavior can be 
traced to the fact that the cavity fields change sign nearly in 
synchronism with the reflections travelling in the coupler. 
Energy coupled from the cavity thus tends to add to the reflected 
wave rather than cancel it. 

To avoid the synchronism between cavity and coupler waves, 
slow wave structures have been tested as coupling elements. The sim- 
plest’ slow wave structure, shown in Figure 5.3, is derived from the 


long strap coupler by zig-zagging it. If this structure is coupled 
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Figure 5.1 Microwave Oven Cavity 
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at the end, its characteristic impedance should be equal to that 
of the input transmission line, while if it is coupled at the 
centre, the characteristic impedance may be doubled. Experiments 
show that the slow wave structure gives improved VSWR, heating 
uniformity and matching between modes. The single stub tuner 
also provides an optimum match over the whole ISM band at one 
particular setting. This indicates that with this structure the 
residual reflections consist of the power that is not coupled 


into the cavity and is reflected from the ends of the coupling 


1 ee 
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ey 


Figure 5.3 Stripline Slow-Wave Coupler 


structure. 
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Dee Reactance Coupling Theory 

(a) Coupling Capacitance Model 

The standard method for modelling and analyzing cavities 
employs a resonant circuit for each mode of the paving A 
capacitor and inductor represent stored electric and magnetic 
energy respectively, while the losses in the walls and dielectric 
may be represented by a resistor in series with the capacitor. 
With a loo) as the coupling element, mutual inductance for each 
mode may be calculated between the loop and the inductor of the 
model. For a coupling element which is small compared to a 
wavelength and couples with the electric fields in the cavity, 


we may use the capacitive coupling model of Figure 5.4. D is 


the capacitance between the cavity and the coupling element. 


feud (> ¢ iS 


Pin L se 


Figure 5.4 Lumped Circuit Model for a Sinale Mode 


The resonant radian frequency of the mode n is given by 


= e 
w Y T/L, (C+D) ee a 


n 


The energy stored in electric fields is given by, 
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94 
Analogously, let the voltage on C be ¢ = <E> /d where d is the 


distance between the wall with the coupling element and the 
opposite wall, as shown in Figure 5.1. <E> is a spatially 
averaged RMS electric field; 
RN at el REY af. 
<E>= eras ' (E a9 ) dV (523) 


where a, b, and d are shown in Figure 5.2. D is given by 
Taichi pep (5.4) 
are 
Where Z is the impedance between the coupling element and ground, 


and ¢, is given by 


Pee Sa (5.5) 
oo ATs ga eb) 


and represents the voltage induced in the coupling element by a 
standing field corresponding to ¢' on C. Then C + D is given by 


solving 


2 


No = 5(CtD)¢' L526) 


of 


The value of L may be determined from the resonant frequency and the 
value of C; the latter calculated from the geometry of the cavity. 
Ignoring cavity losses other than dielectric heating the Q of the 


cavity in a particular mode is given by 


ia) 


where V is the cavity volume abd and ,V is the volume of the load. 


P, the absorbed power is given by 


Dems wie Ra le ay (5.8) 
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R may be obtained from the energy definition of Q: 


Qe = 2w N/P = 1/w,(C+D)R (5.9) 


Losses other than the load are present, for example, the 
metal walls have a finite conductivity and may be deliberately 
coated with a lossy enamel for a base load. The wall losses should 
be represented by a resistor in series with L and the enamel load 
by a resistor in series with R and C. The Q of the model is given 
by 


Q (§.10) 


] 
se 
Qy 


2wN 
where Qn? Qe and Ww are © for the load, enamel (or other protective 


THI 
aie 





load) and the wall losses respectively. For high Q and D«C, the 
dielectric losses may also be represented by a resistance Ro in 


parallel with C given by, 


Ry = <E>¢/d2P = Qu L (5.11) 


The results of the analysis for capacitive coupling are very 


14) for mutual inductance. Most applica- 


similar to Collin's results| 
tors have both capacitive and mutual inductance coupling. 

For multimode applicators, each mode should be represented 
by similar L,C,R equivalent circuits; the required number of these 


being equal to the number of significant modes. 


(b) Results of Model 


Solving for the input admittance of the model , 
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Z 
._ SD(1+sRCt+s LC 
Yan ? Faerie (5.12) 


+sRC+s LC'+s-LCDR) 
Where C' = C+D. Rationalizing, the conductance is 
ole esi ee 722 (5.13) 
(1+p")"-(p/Q)“(1-p "8a ) 
and the susceptance is 


(14p")*-(p/Q)*(1-p*p2a2) 


where wn = I/LC', Q = 1/(u,RC'), p = s/w, ,a = D/(CHD) and g= 1-a 


S is the laplace transform variable and for steady state is equal to 
jw, which means that at resonance p = j. If p ~ j, Q>>] and a<<13 and 


letting | wf» we obtain 


SI) ee (5.15) 
mn £(1-u,2)*4u, 479%) Q?R 
and, 
‘ ~, : 2 

B, (ju) = jwD+jG;Q(u,“-1) (5.16) 

At resonance, these reduce to 
oe 
“in, = “he jwD (5.17) 


Eq. 5.15 shows that if the voltage were kept constant (which is not 
the case for probes on transmission lines) the power absorbed is 
reduced by 3 dB at frequencies where (u2-1) Q=1. At the same 
frequencies the input susceptance is equal to the conductance if D 


is small. On the slope of the function, the power absorbed decreases 
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twice as fast as the admittance. 

The reflection coefficient has a wider bandwidth than the 
input admittance due to the fact that both the incident voltage and 
the reflected voltage are applied to the load. Solving for the 


reflection coefficient, 


fa Gn s JBs,, € G 
$ 
Gin iu JBin Bo 
145Q(u7=1) = (1492 (u2-1)°)R/,2, 


145Q(ue-1) + (140°(u2-1)*)R/a°Z, ad 
where jwD is assumed to be negligible and Lo = 1/G, is the charac- 
teristic impedance of the inlet transmission line. The behavior of 
the reflection coefficient for loop coupling is very similar. 

The voltage transfer function will be used in section 
9.2c. The tank voltage generated by the input voltage is 
: tank 


¢' /¢.. = ———— 
in ] id 
si). ot tank 


sD(sL+s¢LC'R (5.19) 


sD(sL+s?LC'R)+1+SRC+S<LC 


Reducing this to the normalized quantities, and assuming a<<1<<Q, 


eg melt (5.20) 
in 1+p/Q+p* 


At resonance this gives 


_¢' = jaQ (5.21) 
> 
in 


* Spacing between -3dB points and frequencies where G = B respectively. 
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When the other modes are included we may consider all 
the resulting circuits to be in parallel. The total input admit- 


tance is given by 


. au (1450, (u,2-1)/u3) 
in > L, ai aE (5.22) 
n= u, Q, (inte ) JR, 
coitin'!4) uses the same model to analyze the case where 


a loop terminates the transmission line. The result is 


Sa 2 
~~ fo pS L/w i 
£3——=-Gul—+—}. (5.23) 
at 0 n=l w \2,.W 
La a +j /wQ, 
n 
where a, is M/L, and M is the magnetic flux of the n' th mode 


Q 
which Pireade the loop at unit current from the transmission line. 
This expression may be rationalized to give a reflection coefficient 
similar to that obtained for the capacitively coupled cavity. LG is 
the self inductance of the coupling loop and, if it is small, a 


single mode at its resonant frequency gives a resistance of 


ch ae (5:24) 


(c) Distributed Coupling 


The model of Figure 5.4 can be extended to deal with long 
strap couplers and slow wave structures. The coupling capacitance 
and mutual inductance become distributed as shown in Figure 5.5 
and the coupling structure is equivalent to an infinite number of 
such sections of infinitesimal length. Let z be the distance from 
the inlet port of the transmission line measured along the path 


taken by a wave in the coupling structure. Let the coupling 
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capacitance per unit length be given by D k.(z) and the mutual 
inductance by M k,(z). D and M are the averaged magnitudes over 
the coupling structure such that k, (z) and k(z) average to zero 
and bear the z dependence of the coupling. L,C' and R are 
defined in Section 5.2. C and L are, respectively, the self- 
capacitance and self-inductance per unit length of the coupling 


element. The coupling element itself is assumed to be reflec- 


tionless except at its ends. 





Fagure 5.5 Distributed Coupling 


The current iAz through the distributed coupling capaci- 
tance gives rise to two equal voltage waves travelling in the +z 


directions and with magnitudes given by be = > = Az Z/2 where 


is 
Lo is the characteristic impedance of the coupler. Mutual 
inductance generates equal but opposite voltage waves in the 
coupler having a magnitude |wI, az > k,(z)|- The stored magnetic 
energy in the cavity is almost equal to the stored electric energy 
and is given by Me Ger which defines I. Excitation of the cavity 


is modelled by current flow through Dk .(z) and voltage induced in L 


by the coupler current and Mk, (z). The phase reference for the 
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resonant circuit is the phase between its voltage and that at the 


input port of the coupling structure. If the Q of the cavity is 


high, the fields collapse almost simultaneously throughout the 


volume of the cavity. 


Capacitive and inductive coupling generates a voltage in 


the equivalent tank circuit given by 


28 = Ky kz) (o(z)-8') + Kok, (z)4(z)Z, (5 
where ¢(Z) = d¢(z) + o,(z) and i(z) = (4¢(z) - ,.(Z))/Z, 
The constants K, and Ko are given by 
2 
1+p/Qtp 
and 
_ SM ,_1+p/Q 
Be eae. (5 
o 1+p/Qtp 


The coupling also induces waves in the coupling structure given 


he : : 

siose Kk .(z) (4 -~(z)) + Kak, (z)¢ : (5 
aa (z) 9' (5 
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The constants Ky and Ka are given by 
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Equations 5.25, 5.27 and 5.28 form a set of three 
equations in three unknowns which define the voltages resulting 
from a specific excitation of the coupling structure. Integrating 
this set over z and including the excitation and a reflection 
ry of +] at the end, Zp> We obtain: 

z 


z z 
bei Wb ec olen. Pak t 
: L, az 9 So (Ky k +k ky )o_dz F bs (Ky ko-Kok,) ond, 


Zt 
zie K)k.o'dz (5.30) 


vel") = (0) exply(z'-z))+ f° 2% exp y(z-2')), 


az 5 
(5. 31) 
L = 1 zi 2b, ' 
on(Z') = Tebe(z,) exp(y(z ie az exply (2-2 %)) gs 
(oroe) 


where z indicates a generation location, z' a measurement location 
and y is the propagation constant of the coupling element. 
Obviously, these equations are coupled and cannot. be easily solved 
except by difference equation methods. The values of D, M, Ko and 
kK, must be obtained either from a boundary value solution of the 
geometry chosen, or experimentally. 

In Eq. 5.32, the phase angle of all terms changes rapidly 
with z. The integral is much smaller than that of the absolute 
value of the integrand, unless the k's change sign nearly in 
synchronism with the real part of exp(y(z'-z)). In other words, 
phase differences of backward travelling waves in the coupler 


cause cancellation of reflected power unless some of the waves in 


the cavity are nearly synchronized with those in the coupler. 
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For a lossless cavity (R=0), ¢' increases until the unattenuated 
forward wave is reflected back toward the coupler. If R # 0, ¢' is 
lower and the forward wave is attenuated with distance. This may 
be verified by examining the effect of b¢ in Equations 5.30 and 5.31. 
Its phase is synchronous with exp(y(z - z')). 36! will be reduced 


dz 
if »' is large and oe is increased by be and bn: 


Alternatively we may examine the energy transfer. If the 


k's are slowly changing, voltages of nearly equal magnitude and 
phase are induced in a long (but not necessarily straight) section 
of the coupling structure. Because of distributed arrival times at 
the inlet port, the reflected waves cancel and energy does not 
readily escape from the cavity once it has been distributed inside 
the high Q structure. However, there is a gradient in the forward 
wave propagating in the coupling structure and energy can be coupled 
into the fields unless ¢' is large. The energy is randomized and 
cannot be recovered unless the integral of Eq. 5.32 is made non-zero. 
The alternative to slow wave coupling is fast wave coupling, 
where the phase velocity is made very high. This means that the k's 
change sign much more rapidly than cos(y(z-z')) and the reflected 
power integral is small. It is thus possible to use a long straight 
section of waveguide very near cutoff as a coupling structure. It 
would be simply an extension of the transmission structure and would 
have an array of slots to couple power into the cavity. This would 
demand precision construction, but would require fewer parts and 


would take less room than a waveguide slow wave coupling structure. 
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Siegel Results of Slow Wave Coupling 


The basic TEM slow wave coupling structure of Figure 5.3 
has been built and tested. Preliminary measurements of Zo of a 
straight strap shows that a 50 ohm characteristic impedance would 
give little coupling. Several 100 ohm structures tested gave 
much more (but not excessive) coupling in most cavities. The 
presence of capacitance between sections of the structure permitted 
increased spacing to the ground plane at 100 ohms but resulted in 
pass band behavior. In the pass bands, the structure matched well 
to a cavity and to free space; and over half of the 1 to 4 GHz 
range there was less than -6 dB reflected power in free space 
operation. The centre frequencies of the pass bands occur at 
roughly even intervals and are listed in Table 5.1 for several 


experimental structures. The symbols are defined in Figure 5.3. 


Table 5.1 Optimum Coupling Frequencies 


W Ww p f optimum 
10.5 cm Ios -cm Sad bey 2./ GHz ..erce 
8.7 cm ieececm 6.4 cm 2.0 GHz siovurGHzeeecc. 
lore om 1.00 cm 4.0 cm 0.92. GHZ, 295 GHz, etc. 
l26°cK té237cm 6.7 cm 2.45 GHz, etc. 





h = 1 1/2 cm, thickness = 0.2 cm 


Figure 5.6 shows a prototype 80 litre cavity which was used 
for thawing deep frozen, fluorocarbon-DMSO perfused canine kidneys (2°). 
The organ was held in teflon and immersed in 1 litre of FC47 flouro- 
carbon, and placed on a 1" thick nylon turntable. The walls were 


covered with perforated iron sheeting, which increased the losses in 
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Figure 5.6 Experimental Cavity for Thawing Live Organs 





Figure 5.7 Reflections from Experimental Cavity 
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the metal and in combination with the nylon gave a limited Q5° 
This base loading guaranteed good VSWR and limited the maximum 
field concentration in the cavity. 

When the small element and long strap coupling structures 
were used, the reflections from the cavity were unacceptably high. 
With the last slow wave coupling structure of Table 5.1, the maximum 
reflection was reduced from -1dB to -6d and typical reflections 
were about -10dB. The efficiency of heating was about 55% for a 500 
ml water load while for 3 litres of water it was over 85%. The 
system was made useable by the new matching technique; organ heating 
rates of from 2° C to 4°C/sec were achieved in the ranges -79°C to 


55) 


+30°c and +4°C to +302! Uniformity of heating was also improved, 


partly by simultaneously rotating and oscillating the nylon 
disc'*4), 

In Figure 5.7, the reflection coefficient of this cavity is 
shown for a multiplicity of loading conditions as the turntable 
rotates. This and the next two photos were obtained by a frequency 
domain reflectometer and the vertical axis is reflected power, in dB. 
Zero dB means that the energy is totally reflected, and a OdB trace 
was included by removing the cavity completely from the instrument. 
The -10dB line 2 major divisions lower, for example, represents 10% 
reflected power. Similarly, -3dB is 50%, -6dB is 75% and -20dB is 
1% reflected power. The lower a trace is on the photographs, the 
better the cavity is matched. The horizontal axis is the frequency, 
and Figure 5.7 displays the 2.05 to 3.05 GHz band. 


A shield was placed over the coupler in the cavity shown in 
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an attempt to obtain improved matching by reducing the rate of 
coupling to the cavity. As stated, the horizontal and vertical 
scales are 100 MHz and 5 dB per division respectively. The single 
line indicates 0dB, with a marker at 2450 MHz. The shield covered 
40% of the area of the coupling structure and increased the maximum 
reflected power in the 2450 +50 MHz ISM band from -6dB to -4dB. This 
is roughly proportional to the area covered. The load was only 200 
ml and larger loads gave a better match. The bandwidth has been 
narrowed by the single stub tuner but it is still very wide compared 
to the useable bandwidth given by smal] element couplers. The 
failure of the shield to improve the match indicates that the cavity 
is not too heavily coupled, and that the reflections seen are from 
the ends of the coupling structure. 

In Figure 5.8 the slow wave technique and a capacitive 
probe are compared in a 1100 litre empty aluminum cavity. The 
first structure of Table 1 has been used, so the frequency chosen 
was 2.6 GHz +50 MHz. The vertical scale is still 5 dB/division, 
while the horizontal scale is now 10 MHz/division. The capacitive 
probe produced the bottom trace and the slow wave structure the 
upper trace. The 0 dB lines follow the graticule lines with a 
levelling error of about +1 dB. This is a difficult cavity to 
couple to because of its very high Q (and therefore high ¢' 
according to the theory), but there are bands where the slow wave 
structure couples quite well. 

In Figure 5.9 the frequency has been dropped to 1.8 GHz 
where a long strap radiator was optimized. The slow wave structure 


was not quite optimum, but it coupled well. The scale factors are 
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Figure 5.8 Reflections from Slow Wave Coupler and Capacitive Probe 








Figure 5.9 Reflections; Probe, Strap and Slow Wave Coupler 
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the same. From top to bottom are the reflection coefficients 
of the capacitive probe, long strap and slow wave structure 
respectively, in the same large cavity with a 1 litre water load. 
The improvement from the top to bottom trace is typical of the 
behavior of these three structures. 
In this too-large cavity the slow wave structure was 
able to match 11 litres of water to the magnetron with about 
-10 dB average reflected power. The previous cavity was superior 
at the operating frequency and had sufficient base load to match 
well when used with the slow wave coupling structure to heat 
any load. A well designed cavity and the new coupler give a good 
system; if either is poor, the reflections are more severe. 
Striplines are smaller than waveguides, which is physically 
convenient. However, this means that waveguides can have a much 
higher power handling capacity. Several waveguide slow wave 
structures meet the theoretical requirements derived in Section 
5.3. A WR284 slow wave structure was constructed as shown in 
Figure 5.10, with lengths 1 


= 10 1/2 cm, 1, = 3 1/2 cm and the 


1 2 
slots gradually increasing to 0.6 cm in 2 1/2 periods. The free 
space reflected power was less than -20 db of the incident power 
at 2450MHz. When tested in the large 1100 litre cubic cavity it 
gave good matching similar to the stripline structure, although 
the guide wavelength was almost double the free space wavelength. 
The path 1 was large compared to a wavelength. 
It was noted that if the reactance coupling could be 


suppressed, far field effects would give constant radiation 


resistance. With cavities of the size normally found in 
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Figure 5.10 Waveguide Slow Wave Coupler; TE, 


microwave systems, experiments with small elements always 
revealed the mutual reactance (near field) behavior predicted 
in Section 5.2. It is possible to obtain good VSWR and heating 
uniformity in lossy enclosures by using standard high gain 
antennas, but the unused power is not returned to the load, 

and efficiency is usually poor. Other special antennas are 


(54) ,(72) and are not treated here. 


used for concentrated heating 
When the slow wave structure was used to excite a 
single mode in small cavities, it was possible to obtain a good 
match over a wide range of loading conditions. The problem of 
frequency pulling by the load is still present, however, and the 


cavity should either be tuned or be designed for multimode 


operation. This will be treated in the next chapter. 
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5.4 Conclusions and Practical Design Considerations 


The simple small element gives a surprisingly tolerant 
match, and may be used in many situations where the load is quite 
predictable. For example, from Eq. 5.18, a loading range of from 
0.333 to 3 times the optimum R gives a maximum of -6 dB reflected 
power. The 3 dB points for a matched load are about 2.4 times 
as far apart as F,/Q. However, kK. (or k,) varies and the match is 
very dependent on locations of both coupler and load. This is 
overcome by the long strap coupler, giving more useable modes and 
an even wider range of acceptable load variations. The slow wave 
couplers have been shown here to suppress the remaining problems 
of strong reflection for some modes and random matching by the 
stub tuner. 

The slow wave structures used are still rather large, 
approximately a full wavelength in width. The optimization of 
these structures has only been carried as far as selecting geo- 
metries that give pass bands at the desired frequencies. It is 
likely that structures of other dimensions would give further 
improvement in a cavity, although the best match was obtained at 
the frequency corresponding to the optimum free space radiation. 
To check for this, a difference equation computer program for 
Equations 5.30, to 5.32 could be developed. The coupling could be 
increased, since the shielding test indicated that reflections were 
increased somewhat by reducing the effective area of the coupler, 
However, in very high Q cavities, it may be an advantage to shield 


the coupling structure to spread the gradient of the incident wave 
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in the coupling structure over its length. It may be possible to 
parallel more than two sections of coupling structure to obtain 
increased coupling. 

There are no resonances near some frequencies, and 
coupling is poor with any structure. Energy is then reflected 
from the end of the coupling element. It may be that the 
magnetron locks to these frequencies, as indicated by husiae™ 
As he shows, it is desirable that the frequency parameter of the 
load locus of the magnetron Rieke diagram increase while the 
operating frequency of the magnetron decreases. This decreases 
the frequency drift of the magnetron-cavity system as detuning 
occurs. Similarly, by choosing the correct line length from the 
magnetron to the end of the coupling element, the end reflections 
may be placed in a stable operating region on the Rieke diagram, 
or possibly placed so that the magnetron does not operate at this 
frequency. 

From energy considerations, ¢' in a completely lossless 
cavity will rise until total reflection occurs. In low loss 
structures, there is some power absorption, but not enough to 
protect all magnetrons. Then it is necessary to use a permanent 
minimum load. This base load broadens the modes and also reduces 
the reflected power at the centre frequencies. The efficiency of 


Catal 
the cavity is given by Qo/ (9+) where Q.= 0 W and is due to 
e ‘Ww 





the various base loads. The use of the slow wave structure permits 
higher unloaded Q and can thus improve efficiency while permitting 
a system which can be operated unloaded. It should be noted that 


as the Q is raised and coupling improved the fields in the interior 
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are increased and there will be more leakage of energy from the 
cavity unless the door seal is well designed’ *®) , 

A further improvement in heating efficiency has been 
obtained by several workers (32) »(43) by placing the load at the 
focus of an elliptical structure. An alternative is the apex of 
a pyramid. In this type of cavity, ge increases as the focus 
region is approached. The same effect will be obtained in the 


next chapter by adjusting F. 
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CHAPTER SIX 
A SMALL MULTIMODE CAVITY FOR RODLIKE LOADS 


In this chapter it is shown that vanes inserted into a small 
rectangular cavity make multimode operation possible. Tests on 
a preliminary design are given and again the slow wave coupler 


is shown to reduce the effect of load variations. 


6.1 Microwave Heating of Rodlike Loads 


Many presently used microwave threadline driers use a 
single resonant mode which is tuned to the operating frequency of 
the magnetron. This has the advantage that the applicator is small, 
resulting in high efficiency and a wider bandwidth. However, the 
load pulls the resonant frequency downward, and this problem is 
aggravated by small size. A typical cavity is shown in Figure 6.1 (36) | 
In this mode the electric field is parallel to the surface of the 
load and the standard small element coupler is shown. The cavity 
may be tuned by adjusting its volume or shape. In this type of 


system low loss, physically small filaments can be heated, provided 


that the loading does not vary. 
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A load may be considered heavy if it results in a loaded Q 
considerably less than half of the unloaded Q, Q0° Altman (2) 
gives a perturbation analysis for Q and detuning of loaded cavities, 
including cases where the load may not be parallel to the electric 


vectors. The effect of a load on losses is given by 


The subscript o is used for unperturbed variables while the perturbed 
variables are unsubscripted. The denominator is equal to the stored 
energy if the perturbation is small. If the load is parallel to the 


electric field vectors, 


L voltage ot bw (6.2 


fe) é. FE 


The efficiency of heating is given approximately by Q(1/Q-1/Q,) - 
If the volume is decreased, the stored energy decreases and the wall 
losses decrease also, but more slowly. Therefore the efficiency 
increases even though % is reduced. 

This analysis also gives the detuning of the easanthd 


frequency w due to a perturbation: 
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nf Ep is greater than one, as in ordinary dielectrics, the cavity 
resonates at a reduced frequency; while for plasmas, where the 
dielectric constant is less than that of free space, the resonant 
frequency is raised. A small metal object stores little energy at 
high frequencies (equivalent to u+0 and conl?)), and it may be 
shown that a metal perturbation at one of the walls results in an 
increase of resonant frequency if it is located in a region of 
predominantly magnetic energy and a decrease in frequency for 
electric energy. This is the standard method of tuning a cavity 
to compensate for detuning due to the load. The denominator in 
Eq. 6.3 decreases with the cavity volume and detuning becomes more 
severe if the filling factor is increased by using a smaller 
cavity. 


(36) showed that a magnetron can be made to pull into 


Huang 
the operating range of a cavity under favorable conditions. As 
the anode voltage of a magnetron fluctuates due to the power supply 
ripple its resonant frequency changes, often by more than the 
bandwidth of the cavity. If the locus of the load reflection 
coefficient and the magnetron operating point change in opposite 
directions on the Rieke diagram as frequency changes, the system 
becomes more stable. A circulator was used to protect the power 
source, and a reflecting iris and phase shifter to rotate the locus 
of the load. Such a system is expensive, but it may be needed to 
heat a thin, low loss threadline effectively. 

There are several other techniques in use for rodlike 


loads where ey AV is large. A travelling wave applicator can be 


made by generating an electric field component parallel to the 
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axis of a waveguide and passing the filament or rod through it 


axial ly(©5)»(17) This avoids tuning problems but can lower 


efficiency. Double ridged waveguides can be used where the 
load is heavy enough so that the electric field need not be kept 
parallel to the Toad (8) »(41) | Elliptical cavities have been 


designed to control 52(32).(43) | 


Microwave ovens have been used 
for tests but Bhartia et al have correctly noted that these 
cavities are not designed to operate with rodlike loads because 
of their size, They selected a single high order mode and 
tuned the cavity to compensate for loading. Using an alternate 
technique given below, multimode operation in a small volume is 
achieved and a high filling factor can be maintained with reduced 
detuning. 

Plasmas usually are rodlike and are heated in a great 
variety of applicators, as evidenced by a large collection of 
microwave plasma references which was published in the Journal 


(39), 


of Microwave Power Examples are the coaxial cavity of 


(20) 


Fehsenfeld and Broida and the TEI On cavity shown in Figure 


(50) At present 


6.2, or the same structure without an iris 
microwave plasmas are mainly a research tool, although ozone 
manufacture has long been carried out in other types of plasmas. 
With small element coupling, the difference in Q of cavities 
before and after ignition of the plasma is so great that the 
magnetron must be protected by a device such as a circulator. 
The applicator of this chapter or a slow wave ste tineaae. can 


be used to overcome this problem. 
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Figure 6.2 TE) 93 Plasma Cavity 
6.2 Metal Plates in a Rectangular Cavity 


The modes of a rectangular cavity can be made spatially 
degenerate but only two of the three orientations of a mode usually 
give good coupling between the electric field and a rodlike load. 

Of these modes the TE.On modes are particularly promising because 
their electric vectors are always perpendicular to one pair of walls 
and their resonant frequencies are complelely independent of the 
distance between these walls. In addition, if several of these 
structures (or for that matter any set of resonant structures) are 
coupled together the composite structure exhibits new resonant 
frequencies. Thus plates may be introduced into a rectangular 
cavity perpendicular to the electric vectors of the desired pattern, 


resulting in a multimode structure. A loop-coupled smal] multi- 
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mode cavity based on this principle is shown in Figure 6.3 and its 


input impedance versus frequency characteristic is sketched in 


Figure 6.4, 





Figure 6.3 Vaned Multimode Test Cavity 





Figure 6.4 Coupling-induced Multimode Behavior 


In the TET On mode the elementary cavities formed by the 
plates carry a TE 0 waveguide mode which is reflected between end 


(56). Resonance occurs at the frequency where 


walls of the cavity 
these walls are n/2 guide wavelengths apart. In Figure 6.5 three 


such cavities are shown coupled at the ends and the centres. There 
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are two ways in which adjacent cavities can interact, resulting from 
combinations of spatial orientation (or phase) of the same mode in the 
two structures. 

If the modes are in phase, as in cavities 1 and 2, the 
presence of vanes does not materially change the field pattern. If 
the orientations are opposed spatially, as in cavities 2 and 3, the 
fields diffract through the coupling holes and cancel each other. 

End coupling shortens the energy path length resulting in an increased 
resonant frequency. In Figure 6.5 there are three resonant frequen- 
cies corresponding to 0, 1 or 2 phase inversions. Since there 

are two ways to produce a single phase inversion, this mode is 

doubly degenerate, and could give two modes with stagger-tuned 


cavities. 





Figure 6.5 Coupled TE, 9] Cavities 


A rat-race resonator is shown in Figure 6.6. In this structure 

the reflection coefficient for a wave travelling around the ends is 
nearly zero. Resonance occurs when the phase of the fields in 

the cavity is constant. This need not require electric field nulls 


at the ends, and the wave can simply propagate around the loop; 
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hence the name. The opposed orientation in the two cavities which 
are coupled at both ends is a rat-race mode with possible strong 
reflections at the ends of the path. Many rat-race modes are 
permitted in the cavity of Figure 6.3, and resonant frequencies 
lower than the basic TELon resonance occur, provided that the 


TE 0 waveguide mode propagates in the elementary cavities. 





Figure 6.6 Rat-race Resonator 


The problem of coupling to this structure is at least 
as difficult as to larger multimode cavities. Small coupling 
elements give the type of behavior described in Section 5.2 and 
cannot adequately handle detuning problems and a wide loading 
range. When a slow wave structure is employed, this problem can be 
alleviated. There is a physical problem in placing the slow wave 
structure in one of the elementary cavities, since it competes 
for space with the load. However, if the coupling structure is 
placed in the coupling gaps as shown in Figure 6.7 there are 
several frequencies at which the system becomes quite insensitive 


to load changes. 
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Figure 6.7 Vaned Cavity with Slow Wave Coupling Structure 


bes Experimental Cavity 


The cavity shown in Figures 6.3 and 6.7 was constructed 
of brass plate. The inside dimensions of the outer cavity were 
width b = 16 cm, depth d = Fy em and height a = 17 cm. The vanes 
extended the full width of the cavity and were only 135 cm deep, 
leaving coupling gaps 2 cm deep at both the front and back of the 
elementary cavities. Seven vanes each 1.6 mm thick, were placed 
equidistant in the cavity, forming 8 elementary cavities 2 cm in 
height. The vanes and two of the cavity walls had 25 om diameter 
coaxial holes centered By cm from one side of the cavity; and 
3 cm long cut off waveguide chokes were used to prevent leakage. 
The off centre location of the holes permitted good coupling to 
the TEy on or TEs on modes. However, the holes were centred front to 
back and the coupling to the TE, 02 modes was very weak. Loads 
consisting of 7, 14 and 19 mm internal diameter pyrex glass tubes 


filled with water were placed through the holes. To lower the Qo? 


the structure was not soldered. 
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The cavity was coupled by a loop, a probe and a waveguide, 
as well as the third slow wave coupling structure listed in 
Table 5.1. When an electric probe was placed in the bottom cavity, 
many modes were observed, including families of modes corresponding 
to the resonances of the elementary cavities and various alignments 
of the fields. However, these modes coupled best at certain loading 
conditions and, as predicted by Section 5.2, the cavity was not yet 
insensitive to loading. When probe, loop or waveguide coupling was 
used in the front face in one of the coupling gaps, the same 
behavior was observed, with somewhat fewer modes. 

The slow wave structure performed much better, however. 
This structure was oriented both perpendicular to the vanes, as 
shown in Figure 6.7, and parallel to them. The 2 cm spacing between 
vanes was the same as that between the long conductors of this 
coupler, and the middle vane and transmission line were both centred 
in one face. Thus for the parallel orientation these conductors of 
the slow wave eel ing structure were centred on the planes between 
cavities. 

Modes of the elementary cavity and mode density for the 
composite structure are given in Tables 6.1 and 6.2 respectively. 
This technique has increased the number of modes between 2 and 4 Ghz 
from 2 to about 50. Because of rat race modes, resonances started 
at as low as 850 MHz and the cavity design provided three modes near 
915 MHz. The other potentially usable frequencies are listed in 
Table 6.3. There are many modes which do not result in good coupling 
but for several families the cavity could be made, without tuning, 


to accept a heavy dielectric load. 
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Table 6.1] Single Resonator Resonances 
TE, 9] - 1.269 GHz TE5 93 eo le ee GHZ 
TE, 99 - 1.953 GHz TE 399 - 3.292 GHz 
TE59] - 2.06 GHz TE393 - 3.808 GHz 
TEo go =e 00 | GHZ 
Table 6.2 Modes per 100 MHz 
f centre Modes f centre Modes 
2.05 4 2.65 4 
2.15 5 De AS: 7 
enc5 5 e.Bo 6 
rately 4 2.95 7 
2.45 2 3305 7 
ed Ae 5 
Table 6.3 Frequency Bands for Test Cavity 
Empty 
Cavity Effective 
Frequency Coupling Frequency 
GHz Orientation Pulling Comments 
Le parallel - mode separation 
wide 
1.97 perpendicular upward 
2505 parallel none usable 
2.36 perpendicular none usable 
Zeoo parallel down 
2.53,2.65 perpendicular none usable 
2.53%2.65 parallel down 
3.0 parallel none usable 
Jeu) perpendicular spread separation wide 


* Efficiency was low because load was centred in null. 
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In the empty cavity the placement of the slow wave coupling 
structure was designed to excite only the lowest frequency member of 
a family of modes, i.e. the aligned mode*. When the loads were used, 
the other members of the family had a reduced reflection coefficient 
because the fields decayed with distance from the centre of the 
cavity. This gave the effect of pulling the operating frequency 
upward to compensate for the dielectric. Since a distributed coupler 
was employed, the effect of the cavity Q was much less than for 
small element couplers and the cavity gave a base load. Asa 
result there were frequencies at which the magnetron would be 
protected for none or all test loads. 

In the test cavity, however, good performance was difficult to 
achieve and factory alignment on a reflectometer would be needed to 
produce an industrial cavity. As indicated in Table 6.3, the load 
still pulls all the resonant frequencies down and the effective 
centre frequency very often is pulled out of the magnetron's band. 
With large gaps and only seven vanes the mode separation is consid- 
erably more than that for a 80 litre rectangular cavity. 

The efficiency could not be measured directly for the 
trial cavity because the useable bands did not correspond to an 
ISM band where high power was available. The mode density prevents 
accurate Q measurements in the useab’e bands but at other frequencies 
it was possible to isolate modes. The unloaded Q's of the cavity 
were between 270 and 850. A 17 mm pyrex tube filled with water 
gave Q's corresponding to efficiencies from 75% to 93% except for the 


weakly coupled TE,02 modes (35%) where the load was in an electric 


* Tt could be made to excite the other modes by adjusting a 
single stub tuner at the input to the cavity. 
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field null. A 7 mm water load was only slightly less efficient, 
indicating that the large mass of water is reducing E well below 


“a 


ES in Equations 6.1 and 6.3. 


6.4 Application of the Vaned Rectangular Cavity 


The small test cavity gives multimode behaviour similar to 
that shown by the cavities used in microwave ovens. A distributed 
coupling structure reduces sensitivity to loading. The number of 
plates should be increased and the coupling gaps made smaller to 
decrease mode separation. (This would have made interpretation 
of tests difficult.) Other configurations of the slow wave 
coupling element should be tested to improve the number of well 
coupled modes in a band and waveguide slot array coupling may 
make construction easier. Deliberate use and splitting of degen- 
eracies could further help. 

This cavity obviously cannot be made as small as the 
To 19 cylindrical structure, and it has a large internal surface 
area which must carry current. Therefore it is not capable of as 
high an efficiency. However, there is a wide range of loads between 
the filaments described by Huang and the sort of loads that are 
heated in 20 litre or larger microwave oven cavities. The single 
mode cavities suggested for this range are very susceptible to 
detuning, and the previous coupling elements are affected by Q 
changes. There is now no need to use a single mode cavity if the 


load gives Ep EF larger than 1/Q, ina TE, 9] cavity. 
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There are bands for which, in the empty cavity, the 
upper member of the family is coupled best. Such a band could 
be used to ignite and efficiently sustain a plasma, which would 


pull all the modes up in frequency. 
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CHAPTER SEVEN 
Conclusions 


Microwave travelling wave web-heating systems were found 
to be limited by poor control of heating rates in the applicator; 
through both the filling factor F and the field concentration 
factor en It was shown that the TEM mode facilitates the control 
of both factors. Of the possible designs examined, the rectangular 
coaxial structure was found to be best suited for many purposes. 
Detailed analysis of this structure was presented and charts were 
prepared for the design of practical rectangular coaxial appli- 
cators. The computer technique developed may be used for other 
TEM applicators by simply modifying the subroutines that assign 
charge and measurement locations. 

The cutoff frequencies of higher order modes have not 
been calculated for the rectangular TEM structlire os Perturb- 
ation analysis has been used throughout, and the analysis should 
be extended to predict the effect of very heavy Toads (13) » (63) » (67) 
Other TEM structures, such as those suggested in Chapter 2 could be 
used to overcome the problem of energy concentration at corners of 
the conductors. When the web is oriented at an angle to the 
electric field, the attenuation may be greatly reduced, and a 
structure to accomplish this was suggested in Appendix II. 

It was shown that although the physical process is 
different, the inhomogeneity in moisture content in a drying web 
can be conservatively estimated by assuming that it is the same 


as the inhomogeneity in the initial heating profile. Design charts 
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have been prepared to assist in choosing an applicator attenuation, 
and the effect of reflections from the edge of the web have been 
examined. The non-electrical variables such as temperature and 
air flow have been related to the applicator design and their 
effect on homogeneity can be predicted theoretically, although it 
was necessary to assume batch heating, which is not very general. 

Cavity applicators were found to be limited by the 
reflection coefficient of the coupling structure. It has been shown 
that a slow wave coupling structure can reduce reflected power, but 
some base load is still needed to protect the magnetron. Several 
Slow wave coupling structures have been constructed and tested and 
one of them was found to improve the performance of a biomedical 
microwave heating system. 

It has been shown that it is unnecessary to use single 
mode applicators for rodlike materials in many instances where this 
has been done to improve efficiency. The alternative small multi- 
mode cavity has only demonstrated the feasibility of such an appli- 
cator, however, and must be redesigned to produce a practical model. 
Enough detail has been given to guide such a design. 

Three patent disclosures have been filed with Canadian Patents 
and Development Corporation in the interest of the University of Alberta. 
These cover the use of the TEM mode in web applicators, the use of 
slow wave structures as coupling elements for cavities and the small 


multimode cavity for rodlike materials. 
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APPENDIX I 


LIMIT OF VALIDITY OF THE PERTURBATION ANALYSIS 


In the first four chapters of this thesis it was assumed 
that the electric fields existing in the empty applicator would be 
the same as those in the loaded structure. This is an example of 
the perturbation method, which Chatterjee and Chatterjee?) present 
by eigenfunction expansions and Altman 2) by vector products. A 
perturbation term may be added to the solution of a known case, 
which must differ only slightly from the actual situation. The new 
term may then be obtained by using the original solution functions 
with the new boundary conditions. 


The variational methods?) 


is an alternative approximate 
method. Here the values of certain integrals are known to be 
stationary. The calculus of variations is used to find solution 
equations which quarantee this. In electrostatic systems one of 
these stationary integrals is the potential energy, and in this 
appendix the change in the potential eneray will be used to see if 
a perturbation is small. 

In an inhomogeneously loaded waveguide the TE, TM and TEM 
modes do not propagate. For small perturbations new modes exist 
which resemble these modes. Some may be classified as LSE(longitudinal 
section electric), where the electric fields are parallel to the web 
boundaries, and LSM modes where the magnetic fields are parallel to 


(63) 


the web boundaries Mixed modes also exist. A small dielectric 


perturbation will leave the electric fields nearly unchanged in the 


(21) 


LSE modes and the electric flux 0 uncharaed in LSM modes f 
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The energy in a wave in free space travels faster than 
in the dielectric portion of a loaded applicator. As a result, the 
energy is refracted at interfaces. In the dielectric, waves 
approaching the interface at an angle of incidence greater than the 


69) 


critical angle’ will be totally internally reflected, with non- 


propagating fringing fields outside the dielectric decaying exponen- 


tially with distance'2>), 


For lightly loaded rectanaular structures 
Vartanian et a1 667) found that the electric fields of the LSE modes 
were nearly the same as the sinusoidally distributed fields of the 
TE modes of the empty wavequide. As the amount of dielectric increased. 
the distribution in the free space portion of the structure at first 
became linear and then exponential. In other words the perturbing 
element is the web in lightly loaded rectangular wavequides while 
for a heavily enough loaded structure the web carries a dielectric 
mode which is perturbed by the metal waveguide in the. fringing fields. 
If the empty applicator fields predict that a large fraction 
of the power is carried in the web, the perturbation is too peavy'?0) | 
However, this does not imply that the resulting eneray integral is 
equally far from the actual (stationary) value. For LSE modes & is 
higher than predicted, but the perturhed dielectric mode may be quite 


similar and the accuracy may still be good. Thus the error criterion 


that follows is conservative. 


“an “a 


Let . Ee 
(epee Sy 1: alae InsE E dS . (Al.1) 
ear. aes oS, 
se ou ek 


where the electric fields are determined in the empty applicator. The 


electric energy stored in the electric fields in one meter of the web is 
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vi ] ’ a * 
a iis, Be b aE dS (Al.2a) 


The energy stored in the rest of the structure, S-AS, is given by 


] a nk 
Ni thee gf ; 
Ape pte net E.E ds (Al .2b) 


Let the perturbed electric fields be given by kE, where the perturbed 
field distribution differs from the original fields only by a constant. 
If the web is parallel to the original fields, the total energy in the 


perturbed fields is given by 
(Al.3a) 


21 2 oe ! oe 
N =k Je, | [E+E dS+(e-1) f E+ 6 ds | 


total S AS 


Assuming that no energy has been lost in launching the wave, the 


original total energy is given by 


—_— 
> 


Ak 
N = eo ba” ods (Al .3b) 


total #2 “o i 
Then, for a web thin enough so that the electric field is constant 


inside it, we may equate Al.3a and Al.3b and solve for k in terms of 


—. This gives 


- ah (A1.4) 
7 ie oe 5 el 
r 
In the resulting value of k is nearly unity, we may assume 
that the analysis is valid. If for example, the value of k is less 
than 0.707, the perturbation analysis predicts that less than half 
of the power will be transported in the web fields. The actual value 


will, however, be more and the error is large. In the TEy mode in a 
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rectanqular wavequide, Ee = 2. For a web consisting mainly of water, 
F should be less than 0.012. Hence for a WR284 waveguide, this 
analysis is limited to soaked webs up to 0.018" thick. 

This simple test applies if the surface of the web is 
parallel to the unperturbed electric field vectors. For the perpen- 


dicular case, 


k = 2 ee [Al 5) 


Z 
J) +e F( 1 - 1) 


le 


which again should be nearly unity for the perturbation analysis to 


be valid. In either case, if excessive error is predicted, the 


solution can be obtained by other mefhods 
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APPENDIX II 
THE TWIN-LINE APPLICATOR 

In Chapter 2 the circular coaxial line was examined and a 
closely related rectangular structure was analyzed. This has also 
been done with the familiar twin-line, and the results are given in 
this appendix. In contrast to coaxial structures, twin-line structures 
do not give a levelling effect; the wetter spots are heated less. 


In twin-line applicators the attenuation can also be rather low. 
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Figure A2.1 Twin-Line Applicator 
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To provide mixed boundary conditions, the web may be 
placed at an angle in the applicator. This would mean that the 
electric field in the web would be only partially shielded by 
the surface charge and the attenuation would increase with the 


angle of the web. This has been confirmed by experiment. 


Twin-line is similar to single wire waveguides!) in that 
its fringing fields extend far from the conductors but are non- 
radiating like near fields. Most of the energy is carried between 
the conductors, and a lossy dielectric placed in this region will 
experience heating. The boundary conditions for the electric field 
in the web of Figure A2.1 placed along one plane of symmetry are 
not the same as for the webs in the applicators in Chapter 2. The 
electric fields are incident perpendicular to the web and the D 


vector is continuous across the interface. 


Examining Figure A2.1, the power absorbed may be determined 
as a function of web thickness t and spacing 2a. Assuming a line 
charge of magnitude q coulombs/meter, the field E at a distance r is 


given by electrostatic theory; 


Ea a9eer (1) 
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In the absence of the web, the magnitude of the transverse electric 


field at this plane of symmetry is given by 


Eedld =the ha Son Al.2 
| Eaeld =the (2 (A1.2) 


where r = va2 + y? . Adding the web to the system, and assuming 
that it causes only a small perturbation, the normal D = cE is nearly 
unchanged and is continuous across the interface. The power absorbed 


per unit length is given by: 


ee coo) he Rete dV (A2.3) 


dy (A2.4) 


The boundary condition implies 





TT so 2 
= ay Zenuae 
Pueb ~ “Ee, ¢ ah Ge) (a2+y2)? dy 
€ " 
= ule, (4) a (A2.5) 
[S 
The characteristic impedance of the unloaded twin-line is 
given by 
Lo = 1/Cc (A2.6) 


where C, the capacitance per unit length, is 


TE 


0 
Cc = —-_-— (A2.7) 
cosh” (£4) 
and R is the radius of the conductors. For R<<2a, cosh”! (2a/R) 
behaves as 2n(2a/R), and changes in a affect the values of q 


negligibly for constant transmitted power. The power carried by 
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the structure is 


ae ge (A2.8) 


0 
at ¢ volts and the charges are given by 


neal (A2.9) 


Solving for power attenuation at 1 volt between the conductors, 








Pp e 
_ —web _ ra MP ampres el 2 
20 = Pe = WBE, onan (A2.10) 
‘Oh 
oS TE 
fa ro ee 0 -1,2a 
a= (5) wlgeg) (Go) Cosh (Re) 
or, 
iA dated Par il FU =|, 7a 
oh ded rae aad (55) Cosh (e-) (A2.11) 


For constant Ly? this permits attenuation control through control of 
(t"/ay. 

The power drops away rather rapidly as one moves away from 
the applicator. If Pag is the maximum obtainable power absorption 
per square meter, 


é 4 


P E 
hn eerkl..fee a (A2.12) 
Pax Ea e (a’ + yy" 

y 


This function is plotted in Figure A2.2. The power concentration and 
non-radiating geometry permit the use of this structure as a television 


transmission line, but at high power, this unshielded line is unsafe. 
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Figure A2.2 Twin-line Power Concentration Effect 


The simple cylindrical twin-line gives attenuation control 
which is as sensitive to (t'/a) as the rectangular coaxial structure 
is to t/c. It does not have the effect of concentrating energy in 
a region away from the web, and indeed all flux lines must pass 
through the web. (Ga and F are meaningless here where the regions 
are infinite.) To control web field strengths we may use the walls 
of the applicator as shown in Figure A2.3. This has the advantage 
that Lo is more sensitive to be than a. To prevent much change in 
Zo we must move the walls as we adjust a, or else change w to 


compensate for this effect. For w/b large, Lo is proportional to 


Vv b/w 
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Figure Ag.3 Shielded Twin-Line Structure 


A simple twin-line applicator was built as shown in 
Figure A2.4. The characteristic impedance was set to 200 ohms to 
facilitate matching to a 50 ohm measurement system. A tee was 
connected to the input, one conductor and a section of line. This 
line was a phase inverter connecting to the other conductor. At 
frequencies where it was n/2 wavelengths long the desired mode was 
launched. A moving detector measured the amplitude of the electric 
fields in the vicinity of the conductors. Water and ethyl alcohol 
were poured into specially constructed glass tanks 1/8" and 1" wide. 
No detectable attenuation was observed and indeed the water helped 
to balance the line and reduce radiation. 

Equation A2.11 predicts this very low attenuation, and 
shows that the applicator must not be used where c, is high. 
Furthermore, it shows that spots with low are selectively 
heated and these are the already dry areas. Due to the continuity 
of 0, the wet spots dry more slowly. This is the exact opposite to 


the desired levelling effect. 
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Figure A2.4 Test Setup for Twin-line Applicator 
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APPENDIX III 


CORRECTION FOR LINE CHARGE APPROXIMATION 


In Chapter 3, if two strips of charge are adjacent, the 
error in approximating them by line charges may be significant. A 
correction term for Eq. 3.10 can be derived as follows. 

The voltage for a nearby strip of charge of density o 


can be evaluated for a length +c, giving 





2. 2 2 
f= ae fo Cine eed + at eeny2 + aya )Sdz,. (03.1) 
Oo 


OTe 
fe) 


where the voltage is evaluated at y = 0 andc +o. It can be 


shown that 





Figure A3.1 Potential Due to a Strip of Charge 


c 
fAn(yxtd + a)dx = sinh” (3) + c({nic}- 1) 


ta(in|c + Yate | - fnfa}- 1) (A3.2) 
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Then 
¢ = poz {a,2na, + a,an(c +Ma3 + c?) - apina, 
- ayen(c tae © C2) a ry) 


a a 
ee sinh” (4) a Se Sinh” | (—) (A3.3) 


Subtracting the approximate value given by Eq. 3.8 and taking limits, 


we obtain 
a,gn(a,) - ap2n(a,) a, ta 
a ] ] 2 2 ] Zz 
ee A Ma tga + eras (A3.4) 


where Fo is the correction term in Eq. 3.8, i.e. 


, +A/ 2+ ¢2 
6 = pot vim ( Fy + on(S “YEAS )) (A3.5) 


C00 
and r = (a, + a, )/2 and w = (a, ~ a). For a set of strips of equal 
width, the corrections to Eq. 3.8 are given in Table 3.1, where n = 1 
for adjacent strips, 2 for the next further away, etc. The correction 
F,/2ne must be added to Ls, for adjacent strips (but not their images). 


When this was done, Ly was only changed by 0.5%. 
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Table 


A3.1 


] 
0.0452287 


2 
0.0016717 


9 
0.0005149 


13 
0.0002467 


Correction for Strip Width 


‘2 
0.0106180 


6 
0.0011598 


10 
0.0004170 


14 
0.0002127 


3 
0.0046687 


) 
0.0008516 


11 
0.0003 446 


yh) 
0.0001852 


4 
0.0026165 


8 
0.0006518 


12 
0.0002895 


16 
0.0001628 
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